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ABSTRACT 


A computer code has been developed to evaluate the space radi- 
ation environment encountered by geocentric satellites. The 
Short ^rbital £lux integration Program (SOFIP) is a compact 
routine of modular composition, designed mostly with structured 
programming techniques in order to provide not only maximum 
efficiency but also core and time economy and ease of use. The 
program in its simplest form, that is, stripped of all modules, 
produces for a given input trajectory a composite integral orbit 
spectrum of either protons or electrons. Additional features 
such as running printout, exposure index, peaks per orbit, per- 
cent time in electron trapping zones, differential spectrum, 
solar proton fluences, and punched output are available sepa- 
rately or in any combination with the inclusion of the corre- 
sponding (optional) modules. The code is described in detail, 
and the function and usage of the various modules are explained. 
A program listing and sample outputs are attached. 
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SOFIP: 


A SHORT ORBITAL FLUX INTEGRATION PROGRAM 


INTRODUCTION 

The need to predict reliably satellite exposure to trapped particle radi- 
ation was realized soon after the discovery of the terrestrial charged 
particle radiation belts, which coincided with the advent of spaceflight. 

A crude "Orbital F_lux j[ntegration" (OFI) code was developed at Goddard 
Space Fli^t Center as early as 1961. Over the years, a large, sophis- 
ticated, and complex OFI system evolved from these early beginnings 
(UNIFLUX^),a system that processed and analyzed the data into several 
different tabular and graphical presentations. 

However, with the appearance of economical minicomputers and the corres- 
ponding emphasis on small software systems, attributes like "compact", 

"short", "fast", and "versatile" became very important. 

In this paper a ^hort Orbital Rux Integration Program (SOFIP) is being 
presented which, excluding the plotting capability, affords all the options, 
choices, and variations of UNIFLUX, but with substantially reduced core 
requirements and running times. 

Two basic guidelines influenced the creation of SOFIP: structured program- 

ming and modularized organization. These were followed to the greatest 
degree possible or desirable. 

A detailed description of the routine is given in subsequent sections, 
including an analysis of the method employed in the determination of integral, 
difference, and differential fluxes. 

A review of the program organization is given in Figure 1 which depicts the 
structure of the fully implemented (complete) code. Logic flow and decision 
branching is shown in Figure 2. 

The arguments of its input and output vectors (variables and parameters) are 
presented and described in the appropriate sections. Code listing and 
sample outputs are attached. 

SOFIP is written in standard FORTRAN-IV computer language. Card decks are 
available from the National Space Science Data Center in the 029 model IBM 
keypunch format (EBCDIC). The cards are labelled in columns 73-80 as 
SOFIPxxx, where the last three columns (xxx) contain the sequential numbering, 
which is incremented by one. 

A comparison of the time required to compile, linkage edit, and execute 
SOFIP is given in Table 1, The data relate to full length and stripped versions 
of the code, and were obtained for both electron and proton runs, by processing 
720 positions in each case. 
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The approximate amount of storage required by SOFIP is given in Table 2, 
for the various parts of the program, including also the environment models, 
in object form. 

All results were obtained on GSFC's IBM 360/91 and, unless otherwise stated, 
using the FORTRAN IV G compiler. 


THE CURRENT ENVIRONMENT MODELS 

SOFIP is designed to use Vette's standardized models of the terrestrial 
trapped particle environment, as distributed by the National Space Science 
Data Center, Greenbelt, Maryland. New models are periodically being issued 
to replace older versions whenever additional data or information become 
available that permit a significant improvement in the environment des- 
cription, or that indicate a change sufficiently important to warrant such 
an action. All models, both for protons and electrons, represent a static 
environment at a given fixed epoch. However, it was possible to infer from 
the data used in their construction a change in average quiet-time flux 
levels as a function of solar cycle. To date, a continuous temporal 
description of this cycle dependence has not been attempted. Instead, 
separate models were developed for solar minimum and solar maximum con- 
ditions for either species of particles. 

Current at the time of this writing are the following models: 

Solar Max Solar Min 

Protons : AP8-MAX (1970)^ AP8-MIN (1964)2 

Electrons : Inner Zone: AE6 (1980)^ AE5 (1975)'* 

Outer Zone: AEI7-HI (1980)® — - 

— — AEI7-L0 (1980)5 

where the numbers in parentheses indicate the specific fixed epoch (year) 
for which they describe the average environment. 

In regards to the outer zone electron models AEI7-HI and AEI7-L0 it should 
be noted that: 

(a) the version labelled "HI" favors Vampola's fits to the OVl-19 
data, while the version labelled "LO" is representative of all 
the other outer zone data sets presently available at NSSDC. 

(b) these models do not reflect solar cycle conditions and should be 
used indiscriminately for both min and max phases. 

(c) they are interim models which recently replaced the solar min 
and max versions of the older AE4. 
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It should also be noted that the inner zone (^lar max) AE6 does not contain 
any "Starfish" residuals because data now indicate that these electrons are 
no longer present.® 


METHOD 


Integral Flux 

The composite orbit spectrum for integral energies gives the total vehicle 
encountered fluxes, averaged into intensities per second, for 30 discrete 

energy levels: T 

S(>E^.) = cAtEJ^(>E,.) 
m=o 

c=24/T*86400 

where the summation is performed for the entire simulated mission duration 
Tjin hourstand includes all fluxes with energies greater than Ei* At is 
the integration step-size in seconds, J^is the instantaneous integral flux 
obtained fronr^heinodeT^fDr^e i^" energy level* and c is an averaging 
factor. Note that At must have values equivalent to integer minutes (See 
also note on page 7). 

Difference Flux 


The difference flux is calculated from the integral flux S(>Ei) for the 30 
programmed energy levels: 

D(AEi) = S(>Ei) - S(>Ei+i) for i=l, 29 

D(AE3o) = S(>E3o) 

where D is the difference flux in units of particles per square centimeter 
per second per energy interval. It is important to remember that E is not 
constant over either the proton or thp piprtrnn <;ppr tra — 

Differential Flux 

Differential fluxes are only calculated when there exist 10 or more non-zero 
elements in the integral spectrum. That is, fluxes must be defined for at 
least ten descrete energy levels. In that case, the differential fluxes 
are obtained from the composite orbit spectrum by analytic differentiation, 
using the averaged instantaneous values of the total vehicle encountered 
fluxes at the selected energies: 


j(=E^) - 


as(>Ei) 
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for i=l,30 



where j Is the 
results 1n the 
spline fitting 
available, the 


differential flux and E and S are the same as above. The 
program output table represent the derivative of a cubic 
procedure. If there are less than 10 non-zero elements 
program bypasses these calculations. 


APPLICATION AND USE 

SOFIP contains two types of program sections: blocks, which are the 

essential parts of the program and must remain unaltered, and modules, 
which provide additional features and may be removed. For each block and 
module, the following points are discussed: 

1. Function (including output produced, if any), 

2. User input or user action necessary, and 

3. Any restrictions, limitations, or other special considerations. 

There are 11 modules in SOFIP, some of which are paired into packages which 
fulfill a single purpose. These packages and the output they produce are: 

MODULES Output 

ORBIT L- ZONE BREAKDOWN MODULE 

PERCENT TIME MODULE Percent Time Table 

PEAKS PER ORBIT MODULE 

OUTPUT TABLES MODULE 2 Peaks per Orbit Table 

GEOMAGNETIC SHIELDING MODULE 

SOLAR PROTON MODULE Solar Proton Table 


To obtain the desired results from one of these module packages, both 
modules in the particular package must be included in the run. If one 
module is included but the other is not, no data will be outputted from the 
calculation which those modules perform. However, such misuse of the 
module packages will neither cause the run to abend, nor affect the output 
of any other portion of the program. 

BLOCK Q: Initialization 


BLOCK 0 performs the general initialization and preparation of the program. 
It contains all " type" declarations (REAL , INTEGER), and all dimension, 
equivalence, and format statements. 

The user has only to be concerned with one aspect of this initialization: 
the selection of the environment model (s) to be used. Lines 39-44 of 
SOFIP contain C0M*i0N statements for all of the current environment models, 
both proton and electron: 



C COMMON /AP8MAC/DESCR(8), LIST(l) S0FIP039 

C COMMON /AE6MAX/DESCR{8), LIST(l) S0FIP040 

C COMMON /AEI7HI/DESCR7(8), LIST7(1) S0FIP041 

C COMMON /AEI7L0/DESCR7(8), LIST7(1) S0FIP042 

C COMMON /AE5MIN/DESCR(8), LIST(l) S0FIP043 

C COMMON /AP8MIC/DESCR(8), LIST(l) S0FIP044 

For a brief discussion of the models, see the section "The Current Environment 
Models" . 

To select a particular environment model: 

1. Uncomment the COMMON statement which relates to the desired model, 
that is, remove the "C from column one, 

2. Make sure all other COMMON statements are commented out, 

3. Include in the deck to be submitted, the BLOCK DATA subroutine for 
the appropriate model. 

Note that protons require only one model per run, while electrons require 
two, one for the inner zone and one for the outer zone. The two electron 
models are needed regardless of whether the trajectory to be processed 
visits only one of the two electron zones. 

Because SOFIP performs calculations for only one particle species in a 
given run, the model(s) for only one species is needed. In other words, 
any one run will use either one proton model (AP8MAX or AP8MIN) or one 
inner zone electron model CAEGMAX or AE5MIN) and one outer zone eTectron 
model (AEI7HI or AEI7L0). The program cannot and does not check for 
invalid combinations of models, nor for BLOCK DATA subroutines and uncom- 
mented COMMONS statements incorrectly matched. These user errors will 
produce compilation, linkage, or execution errors. 

Note also that there is no provision for changing models during execution. 
Therefore, in a multiple-orbit run, all trajectories are processed for 
the same species and model s(s). 

When SOFIP is run for el ectrons, a diagnostic message may be produced 
during compilation warning that the variable "DESCR7" in line 51 of 
SOFIP has already been dimensioned. Do not change this dimensioning; 
it is necessary when the run is for protons. The warning may be ig- 
nored. 


BLOCK 1: Initialization 

BLOCK 1 performs the initialization of quantities which must be reinitialized 
after each orbit of a multiple-orbit run. In this block, the input variables 
are read and subsequently wri tten out. Table 3 gives the input format and 
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r... ..... .................. 


a brief description of each parameter, two of which require some additional 
conmiefTts : 

NRGYLV is the threshold-energy selector. Its value is an index into the 
ENERGY array. The desired value of NRGYLV is most easily obtained by 
looking at the Composite Orbit Spectrum from a SOFIP run for the correct 
particle species, and counting down to the desired energy level. The 
usual values are: electrons, NRGYLV=5 (0.5 MeV); protons, NRGYLV=4 

(5.0 MeV). 

CUTOFF determines the orbit time at which processing is to be terminated. 

If the end of the orbit tape is reached before orbit time reaches CUTOFF, 
the program will proceed as if CUTOFF had been equal to the time of the 
last point read. This will not cause any errors in the program. 

BLOCK 2: Input 

In BLOCK 2, the trajectory ephemeris tape is read. The program can read a 
tape written in either of two modes, BCD or binary. To read a BCD (EBCDIC, 
or formatted) tape: 


1. The tape must have been written with format 6E18.8. 

2. Comment out lines 133-134, 138-139, and 144. 

3. Uncomment lines 132, 137, and 143. 


To read a binary (unformatted) tape: 

1. The first input variable, PSNTIM (see below) must have been 
written in single precision; the other five elements mast have 
been written in double precision. 

2. Comment out lines 132, 137, and 143. 


3. Uncomment lines 133-134, 138-139, and 144. 


For either input mode, each record of the tape must contain the following 
six variables in the order specified: 


PSNTIM T 
PSNLON 4> 
PSNLAT X 
PSNALT h 
PSNB B 
PSNL L 


Orbit time in decimal hours (must start at 0.0) 

East longitude in decimal degrees 

North latitude in decimal degrees 

Geodetic altitude in kilometers above sea level 

Geomagnetic field magnitude in gauss 

Mcllwain's magnetic shell parameter in earth radii 


The input parameter ISKIP controls the number of records ignored each time 
a new point is called for by the program; only each ISKIPth point on the 
input tape is actually used in performing calculations. 


6 



The positional coordinates of longitude, latitude, and altitude are not 
used in the flux calculations; these calculations are performed with the 
magnetic parameters B and L only. Therefore, it is of no significance 
whether latitude relates to a geocentric or a geodetic reference frame 
(longitude is invariant in the two systems). The altitude, however, is 
used to determine the position of physical perigee in the case of eccen- 
tric trajectories. 

NOTE; Do not use trajectories with stepsizes of less than one minute; 
they will cause the program to abend. Also, the stepsize must 
be constant for any one orbit, because the time integration assumes 
that the increment is not a function of orbit position k 


E Flux,-At= At T Flux,- 
i ^ 1 ' 

BLOCK 3; Calculations 

In BLOCK 3 some preparatory calculations are performed and the fluxes for 
the current position are obtained. No user action is necessary. 


Running Printout Module 


The "Running Printout Module" prints orbit and flux data for each position 
i used in the calculations when the input parameter KPRINT is equal to 1. 
Tf the input parameter KPRINT is not equal to 1, only each KPRINTth point 
is printed. The printed quantities are: 


Orbit time 

Lati tude 

Longi tude 

Altitude 

Field magnitude 

Magnetic shell parameter 

Instantaneous flux 

Time integrated flux 

Orbital flux accumulation 


Ti decimal hours 

Xi decimal degrees 

(()i decimal degrees 

hi kilometers 

Bi gauss 

Li earth radii 

Fi at the position i^ (#/cm^* sec) 

F^At integrated over the interval from i to i+1 

J=EFiAt sum of all fluxes encountered to this point 
i 


The first six quantities are the same as those read from the ephemeris tape. 
For the first position of the orbit, only the positional data are printed. 


If running printout is not desired, delete this module. 


Orbit L-Zone Breakdown Module 


The "Orbit L-Zone Breakdown Module" determines the amount of time spent by 
the trajectory in each of the four zones into which magnetic space can be 
divided on the basis of electron trapping: 

1. Inner zone: outside trapping region (1.0£L<1.1) 
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2. 

Inner zone: 

inside trapping region 

(1.1<L<2.8) 

3. 

Outer zone 


(2.8<L<11.0) 

4. 

External (no 

trapping) 

(L>11.0) 


These data will be used for further calculations in the "Percent Time 
Module". 

Note that the "Orbit L-Zone Breakdown Module" must be used in conjunction 
with the "Percent Time Module". If the percent time information is not 
desired, delete both modules. 

Exposure Index Module 

The "Exposure Index Module" describes, for the selected processing energy, 
the radiation exposure in terms of nine intensity ranges, rising from "zero 
flux" through 10'*-10^, 10^-10^, etc., to "more than 10"^ particles per square 
centimeter per second" in increments of one order of magnitude. The overall 
exposure of the trajectory to each intensity range (in decimal hours) and 
the total number of particles encountered while so exposed are recorded. 

If the exposure index table is not desired, delete this module. 

The exposure index is calculated for particles with E ^ ENERGY (NRGYLV, 
species) where NRGYLV is the threshold energy selector variable (see 
section "BLOCK 2" for further discussion of NRGYLV). 


Peaks -Per- Orb it Module 


The "Peaks-per-Orbit Module" determines: 

1. the instantaneous peak flux per period, in number of parti cl es/cm^* 
sec with energies greater than or equal to the threshold energy 
selected by NRGYLV, 

2. the time (in hours) and position (in h-<|>-X and B-L coordinates) 
at which the peak flux is encountered, and 

3. the total number of particles accumulated per period. 

The "Peaks-per-Orbit Module" must be used in conjunction with the "Out- 
put Tables Module 2". If the peak data is not desired, delete both 
modules. 

Geomagnetic Shielding Module 

This module determines the amount of time the vehicle spends in regions of 
space where L > 5, for later calculations in the "Solar Proton Module". 
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The "Geomagnetic Shielding Module", must be used in conjunction with the 
"Solar Proton Module". If Solar Proton data is not desired, delete both 
modules. 

BLOCK 4: Looping 

BLOCK 4 concludes the trajectory ephemeris read-loop. All blocks and 
modules between Block 2 (trajectory input) and Block 4 are executed for 
each inputted point of the trajectory. No user action is required. 

BLOCK 5: Output Preparation 

In BLOCK 5, the calculations for the composite orbit spectrum are performed. 
No user action is required. 

Percent Time Module 


The "Percent Time Module" takes the information stored in the "Orbit L-Zone 
Breakdown Module." , i.e. , the number of times the vehicle visits each of the 
four zones defined in that module, and calculates the percent of total orbit 
time spent in each; this data is then printed. 

If the percent time table is not desired, delete both this module and the 
"Orbit L-Zone Breakdown Module". 

Differential Spectrum Module 

The "Differential Spectrum Module" calls subroutine DSPCTR, which calculates 
the differential spectrum from the total integral fluxes obtained from the 
environment models. 

If the differential spectrum is desired, include subroutine DSPCTR and this 
module. If the differential spectrum is not desired, delete this module. 

Solar Proton Module 

This module calculates the exposure factor (i.e., the fraction of the orbit 
during which the vehicle is not geomagnetically shielded, but is exposed 
to the interplanetary intensities of energetic solar protons) from the 
value stored in the "Geomagnetic Shielding Module". It then calls the sub- 
routine SOLPRO^ , which calculates probabil istic solar fluences at preselected 
energy levels as a function of mission duration t and confidence level Q. 

There are two elements in the solar proton module which the user may wish to 
alter to meet his specific needs: 

1. Mission duration t: (REALM, variable name in code: T) 

T determi nes the time interval , in non-fractional months, for 
which the solar proton calculations are to be performed. The code 
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is preprogrammed for one year mission duration (t=12). If a different 
length of time is desired, edit card 276 accordingly. The permissible range 
of T values is from 1 to 72 months. 

2. Confidence level Q: (INTEGER*4; variable name in code: IQ) 

Q denotes the level of confidence, in percent, which the user 
wishes to assign to the results; namely, that for the specified 
mission duration the calculated fluences are the smallest values 
that will not be exceeded by actually encountered intensities. 

The preprogranmed confidence level is 90%. If a different value 
of Q is desired, edit card 278 accordingly. Permissible values of 
the variable IQ are integers between 80 and 99, inclusive. 

If the solar proton information is desired: 

1. include this module and the "Geomagnetic Shielding Module", and 

2. include subroutine SOLPRO in the deck to be submitted. 

Otherwise delete the "Solar Proton" and "Geomagnetic Shielding Modules". 
•Qutput-Punch Module 


This module produces a card deck containing some of the calculated results. 
Each card, with exception of the header cards, contains a label in columns 
73-78 and a sequence number in columns 79-80. The label indicates the 
particle species and whether the card contains energies, integral fluxes, 
or differential fluxes. The sequence number will reflect the card's posi- 
tion in the particular section of the deck to which it belongs, e.g., 
integral fluxes, or solar proton energies; sequence numbers range from 1 
to 5 (in the solar proton sections, 1 to 4). 

See Table 4 for a description of the punched output. 

If the "Solar Proton Module" was not included in the run, or if the trajec- 
tory was completely shielded geomagnetically, no cards for solar protons 
are punched. 

Output-Tables Module 1 

When the program is used in its simplest form, this module should be included 
if tabular output is desired. It produces a composite orbit spectrum table, 
containing: 

1. Integral energy levels, in MeV, 

2. Average orbit integrated spectrum, in particles per square centi- 
meter per second , 
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3. Difference flux, in particles per square centimeter per second 
per a£» and 

4. A column labelled "differential flux", which contains only zeroes. 

When the "Differential Spectrum Module" is included and if there are ten 
or more non-zero elements given in the integral spectrum, items (1), (2), 
and (3) remain the same, but (4) now contains average differential flux 
values, in units of particles per square centimeter per second per keV. 

For more information on the calculation of either of these fluxes, see 
the section on "Method!'. 

This module also prints two additional, independent tables for (a) the ex- 
posure index, and (b) the solar proton results. The exposure index table 
presents in the first column the intensity ranges, in the second column 
the total duration of trajectory exposure to each intensity range, and in 
the last column the total number of particles encountered while so exposed. 
The solar proton table lists in the header the mission duration TAU, the 
confidence level Q, the number of anomalously large events NALE predicted 
for TAU and Q, and the geomagnetical ly determined exposure factor used in 
the calculations, and presents two columns containing respectively the 
energy levels in MeV,and the total fluence per square centimeter for each 
energy, 


If the "Solar Proton Module" and/or the "Exposure Index Module" have not 
been included in the run, then the corresponding tables do not appear in 
the print out. 

If no tabular output is required, delete the "Output-Tables Module!". 


Output-Tables Module 2 


This module prints the results calculated in the "Peaks-per-Orbit Module,". 
There are nine columns on this table. Column 1 is an orbit counting 
device, based on the period of the orbit. Column 2 gives the absolute 
instantaneous peak flux encountered during that orbit. Columns 3, 4, and 5 
indicate the spacecraft position in geocentric coordinates at which the peak 
was encountered, while columns 6, 7, and 8 denote respectively the time and 
the magnetic B-L coordinates for this event. Finally, the last column 
indicates the total flux encountered during that particular orbit. It is 
advisable to disregard the last line on this table because many times that 
orbit is incomplete and the fluxes or positions shown do not correspond to 
true peaks. ' 

BLOCK 6: Program Termination 

In Block 6, the program returns to the beginning of BLOCK 1 where it checks 
whether there is another trajectory to be processed. If no other trajectory 
is to be processed, the run terminates. No user action is necessary. 
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Running Times of SOFIP for 720 Input Positions on the IBM 360/91 Computer 



Compile 

Link 

Execute 

Total 

■ 

Protons 

Stripped 

All 

.06/. 03 
.11/. 04 

.00/. 14 
.00/. 18 

.06/. 04 
.08/. 06 

.12/. 21 
.19/. 28 

Electrons 
Stri pped 
All 

j 

.06/. 03 
.11/. 04 

.00/. 15 
.00/. 18 

.06/. 04 
.09/. 06 

.12/. 22 
.20/. 28 

Times are given in c 

lecimal minutes. The first 1 

Figure in each entry is for 


CPU, and the second is for I/O. Values are averages of three runs for each 
type and species. 



Table 2 


Core Requirements for SOFIP 



Core required for 
source program 
(in 1000 bytes) 

S0FIP+ 


Stripped 

11.8 

All Modules 

1 

16.8 

TRARA1+ 

1.5 

TRARA2+ 

2.5 

S0LPR0+ 

2.4 

DSPCTR+ 

8.4 

FORTRAN Library Functions 


Stri pped 

22.7 

All Modules 

24.6 

Environment Models*: For proton runs 

26.1 or 26.8 

For electron runs 

35.6 or 37.3 

See table below for size of individual models 


AP8MAC 

26.1 

AP8MIC 

26.8 

AEI7HI 

17.5 

AEI7L0 

15.8 

AE6MAX 

19.8 

AE5MIN 

^ 

19.8 


■•■Corapiled under FORTRAN 6 
♦Compiled under FORTRAN H 

If SOFIP with all modules, TRARAl, TRARA2, SOLPRO, and DSPCTR are 
compiled under FORTRAN H, the results will be a total of 3.3K less 
than these. 
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TABLE 3 


Input Parameters: Description and Format 


Card 

Columns 

Format 

Name 

Description 

1 

1-12 

3A4 

NAME* 

Any 12-character alphanumeric description of vehicle or orbit 


20-22 

13 

INCL* 

Approximate orbit inclination in degrees 


30-35 

16 

IPR6* 

Approximate orbit perigee in kilometers 


40-45 

16 

lAPG* 

Approximate orbit apogee in kilometers 


50-51 

12 

MODEL* 

Number of field model (from ALLMAG®) used in calculation of 
magnetic parameters B and L in conversion of trajectory ephemeris 


60-68 

F9.6 

PERIOD 

Mathematical period of orbit in decimal hours 


70-76 

F7.2 

BLTIME* 

Epoch for which the coefficients of the field model were evaluated 
for the B-L calculations, in decimal years A.D. 

2 

1-2 

12 

NRGYLV 

Threshold energy selector for running printout, exposure index, 
and peaks 


10-13 

14 

ITAPE* 

Orbit tape identifier 


20-21 

12 

NTABLS 

Number of copies of tables to be produced 

1 

30-35 

F6.2 

CUTOFF 

Orbit time (in decimal hours) at which run is to be terminated 


40-41 

i 

12 

ISKIP 

Program will process only every iSKIPth point on trajectory tape. 
ISKIP=1, all points are processed. 


50-51 

1 

12 

KPRINT 

1 

1 

The "Running-Printout-Module", if included in the run^will print 
only every KPRINT th point of the trajectory points read. Note: 
KPRINT=Q will cause the program to abend. 


Format of input: (3A4,7X,I3 JX,I6,4X,I6,4XJ2,8X,F9.6,1X,F7.2 /I2,7X,I4,6XJ2,8X.F6.2,4XJ2,8X J2) 

*Starred quantities are not required for calculations, but are used only for labelling output. 




TABLE 4 


Punched Output: Description and Format 


Card 

Variable(s) 

Format 

Comments 

1 

NAME(3),INCL,IPRG j 
3A4,1X,I2,7' ,15,'- 

APG,MODEL,BLTIME 

' ,I6,1X,'I(#/CM**2-SEC) D(#/CM**2-SEC-KEV) M0D/TM=' ,11,'/' ,F6.1 

Header card containing vehicle identification*, inclination*, perigee*, 
apogee*, the units of the differential and integral fluxes which are to 
to be punched, and the model number (from ALLMAG*) and epoch for which the 
B and L of the trajectory tape were calculated^. 

2-6 

ENERGY ( 30, ITYPE) 
1P6E12.4 

The 30 integral threshold energies, in MeV, for the particle species 
considered in this run. 

m 

AIFLXS(30\ 

1P6E12.4 

The orbit integrated, integral fluxes for the 30 energy levels punched 
previously, averaged into units of parti cl es/cm^’ sec. 

12-16 

0IFSPC(30) 

1P6E12.4 

The differential fluxes obtained from the preceeding integral spectrum 
in units of particles/cro^’sec*keV. Zeros will be punched if the "Differential 
Spectrum Module" was omitted. 

17 ^ 

IQ,T, INALE, EXPFCT 

'SOLAR PROTONS #ENERGIES=20 Q=',I2,' TAU=',F4.1,' NALE=',I1,' EXPFCTR=* ,F5.2 

Header card giving the number of Solar Proton energy levels processed, 
and containing confidence level Q, mission duration tau, number of AL 
events for given Q and tau, and geomagnetic shielding effect in the form 
of an exposure factor. 

18-21 

SPNRG(20) 

1P5E12.4 

The 20 integral threshold energies, in MeV, for the solar flare protons. 

22-25 

F(20) 

1P5E12.4 

The total unattenuated interplanetary solar flare proton fluences for the 
preceeding integral spectrum, modified by the exposure factor punched above, 
and for the Q and tau given in the solar proton header card, in units of 
parti cl es/cm^. 


♦taken directly from input 
















BLOCK 0 

INITIALiZATION OF PROGRAM 38-88 


BLOCK 1 

INITIALIZATION OF DATA 
PARAMETER INPUT 

89-130 

BLOCK 2 

INPUT POSITIONAL DATA 
FROM TAPE (TRAJECTORY) 

131-146 

BLOCKS 
OBTAIN FLUXES 
PERFORM CALCULATIONS 

147-174 

RUNNING PRINTOUT MODULE 

PRINT POSITIONAL DATA 

175-192 


BLOCKS 

OUTPUT PREPARATION 244-254 


CALCULATE PERCENTAGE OF 


TOTAL TIME SPENT IN EACH 


ZONE 

255-271 

DIFFERENTIAL SPECTRUM MODULE 

CALL SUBROUTINE DSPCTR 

272-273 


CALL SUBROUTINE SO LPRO 


COUNT TIMES EACH ZONE 
IS VISITED 

193-201 

EXPOSURE INDEX MODULE 


STORE FLUXES AND TIMES 


BY INTENSITY RANGES 

202-213 

PEAKS PER ORDIT MODULE 


DETERMINE PEAKS AND 


POSITION, PLUS TOTAL, PER 


PERIOD 

214-237 


PUNCH INTEGRAL ANO 
OIFFERENTIAL SPECTRA, 
SOLAR PROTONS 




Wmw 


288-307 


WITH ORBIT INTEGRATEO SPECTRUM 
ANO WITH OR WITHOUT OIFFERENTIAL 
SPECTRUM AND/OR EXPOSURE INDEX 
AND/OR SOLAR PROTONS 

308-381 


tillJUl 


PEAKS AND POSITION PLUS 
TOTAL PER PERIOD 


382-396 


DETERMINE WHETHER 


POSITION IS GEOMAGNETIC- 


ALLY SHIELDED (STEP COUNT 


FOR L<5) 

238-240 

BLOCK 4 

END OF READ LOOP 

241-243 















FIGURE 2: FLOW DIAGRAM FOR SOFIP 



I BLOCK 0 


END OF INPUT - ALL TRAJECTORIES PROCESSED 
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5}{ s;< 5}C 5j? Jjc 3jc ;|c 5{C :<c ^ :{c 5;s ^ ;i{c s;c 5|i 5j! sj: >;< i|c ^ s;s s|« i;c s|c ^ ^ ^ ^ sic 5^ 5|C 5jc V ^ ^ ip S|C 3)5 sp 5,’^ S F T P () 0 1 

#Jic##«>}:#»ic#>is*# SHORT ORBITAL FLUX INTBORATIOM PROGRAK >;! SO F I POO/:’ 

POR USE WITH NS.SDCS STANDARD EMV I R OM^ENT ''1005=1. S 5|/>;«s!?5;j>;csis#SnF f ROOi 
:ic :jc # sl^ic >|c >;<:{/ a:c>;c # >ic :jc ^ sj: ^ >|t jjc >{/ 9|<:{: s|c sic ijc 9|2 :[c :;o|: ^ 9|: 3'^ >)c >|c V ^ >i‘ >i< >i< >r ^ >1‘ >;< >'ic S 0 F I P 0 0 A- 

#* designed and tested by STASSIMGPOULOSt Hpi^ERT, BUTLER* f> kaRTh »i«>;‘SnF T POO b 

** CODE 601, NASA/GODDARD SPACE FLIGHT CENTER: GRFEMBELT, MD. ?077l >l«>?SnF T P006 
*=<= SINGLE PRECISION DECK FOR FORTRAN IV ( EBCDIC, U?9 PUNCH) >:<>i‘SnF I ^00 7 

** TRAJECTORY INPUT FROM UNFORMATTED BINARY OR BCD FORMATTED TAPE ^^i^SOP T POOR 

^ sp s;c sp ^spsjcsp spspsp Jp 5p 3}c3p 5p5p spsp 3p ^ 3p :^cs;t::^p :^c sp 3lc sjc 3;c yp ip 3|csp5p a.t ijcsp Jic 5P ::p5p s;c yp 5p s;c ip ip sp sjc ip pc i{c i;c ;p ^ ^ pc sp sp i;c ip ip ip ip sp sp ^c ;p S Q p J P () f ) 9 


>p3pip 

spspip 


INPUT PARAMETERS: 


SOFIPOlO 
SOF IPOl 1 


ipspsp 

5p 

NAME : 

12-CHARACTER 

MISSION (OR PROJECT) NAME 

SOFIPOI? 

spspip 

a|« 

INCL 

APPROXIMATE 

INCLINATION OF ORBIT PLANE IN D'^GREf^S 

( T>!<4)SnFIP01i 

spipip 

3i= 

IPRG : 

APPROXIMATE 

PERIGEE ALTITUDE U; KILOMETERS 

( T3!/4 ) SOF T P014 

^spjp 


I APG 

APPROXIMATE 

APOGEE altitude IN KILOMETERS 

( T>i=4)SnFlP0l3 

spipip 


MODEL 

NUMBER OF FIELD-MODEL USED IN B/L CALCULATTPM 

(Hv4).SnFIPf.U6 

SpipSp 


PERIOD 

MATHEMATICAL 

PERIOD OF ORBIT IN HOURS 

( R’;=4)SPFIP017 

spspsp 

5p 

BLTIMEJ 

EPOCH OF FIE 

LD-MODEL USED. IN B/L CALCULATION 

(P=X4).S0FTP018 

spile sp 

5|c 

NRGYLV 

THRESHOLD-ENERGY SELECTOR FOR RUNNING PRIMTOUT 

( ) SPF I P019 

3p5p5{c 

❖ 

ITAPE 

8/L ORBIT TAPE inENTIFIER, < IDOOO 

( I=i=4)SUFTP0-20 

spspsp 


NTABLS 

5 OF OUTPUT- 

TABLE- SETS PER TRAJPCTORY 

( T>:‘4)SnFIP021 

spspsp 

3i« 

CUTOFFS 

ORBIT DURATION IN DECIMAL HOURS 

( R’l-A ) SUFI P022 

spspsp 


I SKIP 

POSITION SKIPPING CONTROL 

( T--i54 ) SOF IPO? 3 

5!«5;<3p 


KPRINTs 

RUNNING PRIN 

TOUT CONTROL 

{ 1*4 ) SOF I P024 

s|c^c^ 





SOF I PO?b 

^ 

INPUT VARIABLES: 


SUF I D0?6 

spile sp 

'I' 

PSNTIM 

1 POSITIONAL 

TIME ( DECIMAL HOURS ) 

S0FTP027 

i{c sp sp 

:5c 

PSNLOMs 

n 

LONGITUDE (DEGREES) 

SUF T P028 

vV 

‘V* 


PSNLAT 

• ii 

LATITUDE ( degrees ) 

^>CiF I P02 P 

spspsp 


PSNALTs 

If 

ALTITUDE (KILOM-TERR) 

Sf)F I P0 30 

.spipip 

❖ 

PSNB 

• n 

FIELD magnitude ( GAUSS ) 

SUF I P031 

spspip 


PSNI, 

11 

SHELL parameter (eaRTM RADII) 

SOFT P0 32 


3itJ!cj!s;:s3;!5}::{!#>j::{c3;t3!s3|s5;c:;; 


' 3!/ 3jc 5{c -ic sic 3|? 5|s sjc ;Jc J|s i[c 5|s 3S sji i[s ;^c s:/ 5{s 3J; ;;c 5;/ 3 


C 3{C 3iC 3jC 3|C 3|C 3|C SjC SjC SjC sjc 3|C SjC 3|C 3jC S|C 3|C 3JC Sjc 3|C 3*^ 3)C 3|C 3,/ S t 1 I P 0 3 3 


*=!=* TO HEAD BCD FORMATTED ORBIT TAPES, UNCOMMENT LIMES 13?, 137, & 1^3, 
**3ic COMMENT OUT LINES 1 33- 1 3 A , 1 38-1 39 , f, 144. 

**3ic TO READ UNFORMATTED BINARY ORBIT TAPES, UNCOMMENT LINFS 133-13^, 
*** 138-139, & 144. COMMENT OUT LINES 13?, 137, S 143, 


3{c:{C3iC 3!c3ic^3:c3!c3;c3ics;c3ic3lc3;c3;csjc3;c3 


3|C 3}C 3^ 


BLOCK 0: INITI ALT 7ATI ONi 


ip ip ip ip s,c 




(V ,y, ,f, 


COMMON /AP8MAC/DESCR (8) .LIST! 1 ) 
COMMON /AE6MAX/DESCR(8) ,LIST( 1) 
COMMON /AEI7HI/0ESCR7(8), LIST7I 1 ) 
COMMON /AEI7LO/DESCR7I 8) ,LIST7( 1) 
COMMON /AE5MI N/OESCR ( 8 ) , LIST! 1 ) 


SOF T P034 
S n F T PO 3 b 
S0FIP036 
SOF I P037 
=:=SOF I P038 
SOF I P039 
SOF I P040 
SOF I P041 
SOF I P04? 
SOF I P043 



CnKMflM /aPBMIC/OESCrt ( H ) , LI ST( 1 ) S0FIP044 

(4. 7)/ 'HcMORICK ' , ' SKCAIK! 9't'9-TERH G't'SFC 9/B5'.' CSnFTP045 
;*>AIN FT't'.AL. 1?0','-TEKM GS','FC CA I NSL A ' . ' NGE L 143 ' . ' - TSOF I P046 

S.ERW P()','(',0 10/48'. '• CAINaSW , ' EENEY 1 ? ' . ' n-JERM P'.'OGn 8/49',' S0FIP047 
4 IGRF',* 1985.:) '.'80-TERH ','10/48 ',' I.FATONi ','MALIN EV ' , ' AMSOF I P04B 

<bS 80-T','ERM 1945',' HIJkKI'.'TZ US CE','GS 14fi-T','FRN 1970'/, S0FIP049 
s<rtnOLrtL*8 (4 ) , APB/ ' APB ' /,MAX/'MAX '/,NIL'/'M1U '/,LnW/'S LO'/, SOFIPOSO 

SHAG/ ' '■lAC ' / ,DESCR7( H ) . B I UDMY»8 ( 5 ) , AnUN’Y.Y ( 4 ) , M0n7/' LO-7' / SOF IP051 

Ri^AL E"ERGY( 31 ,2 ) /2 . , 3 . , 4 . , 5 . , 4 . , H . , 1 0. , 1 5 . , ? 0. , ? 5 . , 30 . , 35 . , 40 . , SOF IPO 52 
4A.4. ,30. , 5 5. ,40. , 70. ,80. ,90. , 100. , 125. , 150. , 175. , 200 . , ? 50 . , 300 . , SnFIP05 3 
4350. ,400. , 500. , 0. , . 1 . .2, .3, .<i, .5, .4, .7, .8, .9, 1 . , 1 .25, 1. 5, 1 .75,2. , S0FIP054 
S2. 25, 2. 5, 2. 7 5, 3.. 3. 25. 3. 5, 3. 75, 4., 4. 25, ^.5, 4. 75, 5., 5. 5, 4., 4. <^.7., SGFIP0 5 5 
40.0/, SPURG( 20) /in. , 20. .30. ,40. , 50 . , 40 . , 70 . , 80 . , 90 . , 100. , 110. , 120 . .SOF IPO 56 
4130. . 140. , 150. , 140. , 1 70. . 180. ,190. ,200. / S0FIP057 

INTEGER MRGRi;G( 1 0 , 2 ) / 1 , 3 , 5 , 7 , 1 2 , 20 , 2 2 , 24 , 30 , 3 1 , 1 , 5, 8 , 10 , 1 2 , 1 3 , 1 ^ , SOF IPO 58 

4 22,30,31/, I/nM= ( 120) /Ki*\, 17*2,93x<3/,''iRBITn/l / S0FTP059 

n I MENS KIN FLUXES ( 30) , ALGFLX ( 30 ) , ALNF L X ( 30 ) , 0 I FS PC ( 30 ) , F XPF L X ( 1 0 ) , SOFT PO 60 

4PKyALU( 50,8 ) , AI FLXS ( 30 ) , EURUGS ( 1 1 ) , EXPTI M( 10 ) , I YPn( 3 ) , LCfUIMK 4) , SPFIP041 
^niEFi.X ( 30) AYE ( 3 ) ,PKFLX( 50) ,PKTIRU 50) , PTl MF( 4) , S0FIP042 

•kPKLUiv ( 50 ) ,PKLAT( 50 ) , PKALK 50 ) , PK8( 50 ) , PKL( 50) , TAUFLXI 50) , F( 20 ) SOF IP04 3 

EOIJIVALEMCE ( PK V ALU ( 1 , 1 ) , PKF LX ( 1 ) ) , ( PK Vfi I.U ( 1 , 2 ) , PKLON( 1) ) , ( PRVALU ( 1 SOF I PQ64 

5 ,3 ) ,PKLAT <!)),( PKVALUI 1,4), PKALT( 1 ) ) , ( PKVAIJ.H 1,5), PKTIMC 1 ) ) , ( PKVALSQF I P045 

SIH 1 ,6) ,PKH( 1 ) ) , ( PKVALLU 1,7) ,PKL( 1 ) ) , ( PXyAl.Ul 1,8) ,TAUFI.X( 1 ) ) SOFT P046 

real TYPLRL ( 3 , 2 ) / ' PR ' , '0701'" , ' S ' , ' FLFC ' , 'TROU ' , ' S MI'/, SnFTP047 

4F IKMGS ( 11 ) / ' ■). eo' , ' 1 .E0» , ' l.El • , • 1 . E2 ' , ' 1 . E3' , • 1 . E4' , ' 1 .F5 ' , SOFT P048 

S*1.E4','1.E7',' nVER ' , ' ' / , XLABFL=i=8( 3,2 ) / ' PRNRGY' , ' PR lUTG' , ' PRO I F F ' SOF I P049 


■b, ' ELNRGY ' , ' ELIOTG ' , ' FLOI FF ' / ,PRnTL»^ 8 ( 2 ) / ' SPNRGY' , ' SPFLI!X' / SOF I PO 70 

3 Fi'.HMAT ( ' 1HAH6 = '.3A4/' 1UCL = ',13/' IPRG = ', 14/' lAPG = ', S0FIP071 

SI4/' ITAPE = ',14/' MontL = ',12 /' PEkIOD= ',F9.4/ SQFIP072 

4' RLTIME= ',F7.2/' NRGI_EV= ',12/' NTABLS= ',12/' CUTOFF= ',F4.2/ 40810073 

S' ISKIP= ',12/' KPRINT= ',12//) SnFTP074 

2 FOkMAK '] ' ,131( '=X' )/' # SOFIP : SHORT ORBITAL FLUX INTFGR. PROGR A^SOF I P07 5 

S FOK STANOARO ■■•ssnc PRofniv and electron EMVIR. MODELS (SPFCIFS CnUSOFl P074 
SSIDEREn SEPARATELY) *'/' =? MAGMETIC PARAMETERS B AND L COMPUTED MTS0FIP077 
STH GEOi-'AGN. FIELD MODEL', 13,': ',4AH,' =? COFFF. UPDATED TO : ' , F7 . 1 , SOF I P07B 

S' # PROJECT :',3A4,i INCLIN= ' , I 3, ' OEG ^ ppr J 1 -;= ' , i 5 , ' KM APSGFIP079 

S0r;= ' , I 4, 'KM =!= B/L TAPE = TD' , 14, ' P ER I 0D= ' , F7 . 3 , ' HRS ^ SOLAR ' , A3 , SOF T P080 
S'lMUM *'/' * FOK INFORMATIOM OR FXPLAMA T TOW COWTACT F.G. STASSSOF I P08 t 

SIWOPOULUS AT MASA-GSFC,CODE 40T, GRFEN 8 FLT, MARYLAMD 20 7-71 , TEL . ( 3SQF I P082 
SOI )-344-8067 =!='/ 1 X, 13 I ('=:?')// ) S0FTP083 

3 FORMAT ( 2 ( ' I ' /12 ( ' 0 ' / ) /53X ,28{ » =X' ) / 53X , '=!^=i= ' , 4X, 3A4, 4X ,' '/ 5 3X S.OF I P084 

S, 13, 'DEG/', 15, 'KM/', 14, 'KM ** ' / 5 3X , 28 ( ' ' ) / ) ) , SOFIPOBS 

-Y FORMAT ( 3A4,7X, 13, 7X, I4,4X, I4,4X, I 2 , 8 X, F9. 4, IX, F7.2/ T 2, 7X, T4, SOFIPO 86 

s4X,I2,8X,F4.2,4X,I2,8X,I2) SOFIP087 


5 FORM AT(4E18.H) 


SOF 1 P088 


L-2 



o o 


{j ^ 'I' 'fi 
] > 


1 S 
17 

20 


( MAi'^E. INCL, IPRG, I APG, 1 = 1 , ? ) 


RI.OCK 1: I W I T I A L I Z A T I DN 

Rh AD ( 5 t-A, ENQ=909 ) IV AME , J 9CL , I PRG , I APG , MODEL » P ER I 00 ♦ BLT I MF » NRGYLV/ , 
S' I TAPE , NT A BLS , CUTOFF , I SK I P , KPR IMT 
1 TYPF=1 
9 hITF(6,3) 

MmBIT = l 
IPASS=1 

ipkimt=<primt 

ASSIGN 110 TO ?JGn2 
L = 0 

LSIJN = 0 
EXPFCT=0.0 
XAMivII-' = MAX 
ISWTCH=1 
iFioESCKi 1 ) .eo. 

ITYPF = 2 

assign 120 TO NG02 
IF (OF SCR (2 ) .NE.HAX) 

IP(0ESCR7(2).Ei 
GO TO 17 

IF(DESCR(2).NE.MAC) 
no 20 1=1,4 
LCOUNT ( I ) =0 
MnOLRL(I) = moolabci 
TAU = PEKlOn 
FLXSO’'' = 0.0 
OF LX SO = 0.0 
PPAK = -1.0 
DO 3') NRMG=1,10 

EMRNGS(NRNG) = ENERGY ( NH GRNG ( NRMG , ITYPF) , ITYPE) 

EXPTIM(NRMG) = 0.0 
30 EXPFLX(MRMG) = 0.0 


I PH) GO TO 15 

XAHNIM=MIM 
!.M0D7) TYPLRL(3,2)=L0W 

X f’'* M I f'*' = tA I M 
MOPFL ) 


■^SOF I P089 
S0FIP090 
S0FIP091 
S0FIP092 
SOF IP093 
SOFT PU94 
SOF IP095 
SOFT P096 
S0FIP097 
S0FIP098 
S0FIP099 
SOF I Pi 00 
SOF IPlOl 
SOF I Pi 02 
SOFIP103 
S0FIP104 
S0FIP105 
S0FIP106 
SOF I P107 
SOFI P108 
S0FIP109 
SOFTPllO 
SOFIPlll 
S0FIP112 
SOFIP113 
S0FIP114 
S0FIP115 
S0FIP116 
S0FIP117 
SOFIPl 18; 
S0FIP119 
S0FIP120 


on 35 IMRGSP = 1,20 

S0FIP121 

35 F (NKUSP) =0.0 

SOFIP122 

DO 40 NRG=1,30 

S0FIP123 

AIFLXS (NRG) =0.D 

S0FIP124 

ALNFLX ( NRG) = 0.0 

S0F1P125 

niFSPC(NRG) = 0.-'^ 

S0FTP126 

^0 FLUXFS(MRG) = o.O 

S0FIP127 


WRITE OUT INPUT P AR Afv-ETER S S0FIP128 

NR I TE ( H, 1 )MAMF, I NCL , I PRG , I APG , I T AP E , MODE L i PER I 00, BLT IME, NRGYLV, SOF I Pi 2 9 

■!>NTAHLS .CUTOFF, I SK I P, KPR INT SOFI PI 30 

BLOCK ?; IMPI IT I PI 31 

RE AD ( 9, 5 , EN0=4 )0, ERR = 10) PSNTMl , PSNLNl , PSNLTl , PSNAL 1 , pSNBI , PSNL l SOFI PI 32 

R E A 0 ( 9 , E N D = 4 0 0 , E R R = 1 0 ) P S I'l T M I . P S N 1. 1' 1 , D U M M Y , P S N L T 1 , D U M M Y , P S M A L 1 , S 0 F I P 1 3 3 


L-3 



5>DUMMY , PSM8 I , OU '-''MY , P SNL 1 , nUMKY 
TMLAST = PSNTKl 
80 on 60 ISKP=1 , ISKIP 

C RE AD { 9, 5, END = 400, ERR=10 ) P SNT I 0 , P SML 00 , PSi'il. aj, PSOAI..T. PS'Mrt , pSmi. 

REAO( 9 , END=4C0, ERR = 10) P SOT I M , P SNLOW , DUMOY , P SNO AT , 01 !MM Y , P S^''- LT , 

S DU M Ni Y , P S N 8 , 0 U H M Y , P S lY L , D U M O Y 
60 CGNTIMUE 

IF(PSNTIM.LE. CUTOFF) CO TO 65’ 

C DUMMY READ LOOP TO READ TO END OF FILE 

C 66 READ ( 9, 5, END=400, ERR=10 ) ADUHPY 

66 READ(9,END=400,FRR = 10) 8 I N T I M , 8 I NDM Y 

GO TO 66 
65 CONTIRiUE 

C ^ ^ V ^ ^ ^ V 3(S SjC 5gC ^ 3|£ 5(C sj! l[> 5|I sjs l[i Sjt 6 L 0 C R 8 • C A I C U L A T I i ! 9 ^ V ^ ^ ^ 'P ^ ■ 

c CALCULATE KPSTEP (NUMBER OF MINUTES BETWFF'v POINTS ON B/L T^^P^) 

GO TO (70,80), I PASS 

70 KPSTEP = INT (( PSNT IM-TMI_aST) /. 0166667 + 0 . ) ) 

80 TMLAST = PSNTIM 

C TF-ST L-VALUE C BYPASS FLUX CALCULATIOMS IF WARRAMTPD 

IF(PSNL.GT.O.O.AND.PSML.LT. l?-0) GO TO MGO''' , ( M 0, I ?0 ) 

DO 100 NRG=1,30 
100 FLUXES! NRG) = 0.0 
GO TO 170 

C *** OBTAIN COI^MON LOGARITHM OF POSITIONAL FLUXES (ALCFI.X) 

C PROTONS 

110 CALL TRARAl ( D ESCR , L I ST , P Sf'M.. . P SN B , F MFR GY ( 1 , 1 ) , ALCFLXI I ) ,30) 

GO TO 140 
C ELECTRONS 

120 IF( INT( l()0.0-'PSNL + 0.2) .LE.280) GO TO 130 

CALL TRARAl ( DESCR7, L I ST7 , PSNL , PSNB, Ei^'FRGY ( i , 2 ) , A L GFL X ( I) , 30 ) 

GO TO 140 

130 CALL TR ARAl ( OESCK , L I ST , P SN L , P SNB , FW ER GY ( I , 2 ) , A L GFL X ( 1 ) , 30 ) 

C CONVERT LOG-FLUX TO FLUX 

lAO DO 150 NRG=1,3) 

FLUXES (NRG) = 10. Of- ALGFLX ( NRG ) 

150 IF(FLUXES(MRG).LT.l.noi) FLUXES(NRG) = O.d 
C SUM FLUXES FOR (A) RUNNING PRINTOUT, (R) TABULAR OUTPUT 

FLXSUM = FLXSUM + FLUXES(NRGYLV)'+FLOAT ( KPSTEP ) *60, 

DO 160 MRG=1,30 

160 AIFLXS(MRG) = A I FL XS ( NRG ) +F LUXE S ( NR G ) 

170 CONTINUE 

Q :(!>;: RUNNING P R I iV TO U T I'-'n DU I E 

GO TO (200,210), I PASS 

200 WRITE (6,2) MODE L , MODI BL . BL T I ME , N AME , INCL, IPRG, lAPO, ITA pf, PERIOD 
SXAMMIM 


S f !• r T 1 3 4 
SOP T Pi 3 d 
SOFT iJl 3f- 
SUE f Pi 37 
SHF IFl 38 
SPF 1 Pi 39 
SOF I Pi 40 
Sn- I P141 
S n F I H 1 92 
SOF I PI 43 
SnFTPl99 
S 0 F T P 1 4 5 
SOFT PI 46 
>?SOF TP147 
SOF I P 1 48 
SOF I Pi 49 
SOF I Pi 50 
S n F I P I 51 
SOF T Pi 52 
Slip T PI 53 
SUE I PI 54 
S n F I P i 5 5 
SOF I Pi 56 
SOF rPl 57 
SHpi Ml 5« 
SOF I Pi 5 p 
SOF 1 Pi 60 
SOFTdi A 1 
SOF I Pi 62 
SOF I P163 
SOF I Pi 64 
S0FTP165 
SOF T Pi 66 
S0FTP167 
SOFI P168 
SOF I Pi 69 
SOF I Pi 70 
SnFTPl71 
SOFT Pi 72 
SOF I Pi 73 
90F I PI 74 
=i=SnFIPl7 5 
S0FTP176 
SOFIP177 
SOF T Pi 78 
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c 

c 

c 


c 

c 


c 

c 


WRITE ( 6, 201 ) (TYPLBLT I . ITYPE ) , T = 1 , 3 ) , EM ER GY ( NR 0 Y|. V , T T Yp P ) -SOFlPlTW 

201 FORMAT! 'O' ,21X, , 3 A4 , ' ( E> ' , G9 . 3 , ' ME V ) // i LnMr,snF 1 Pi «0 

S. LAT, ALT. FIELO LINE ORBIT POSTTinWAL T T (■■iE- U T PSivF I P 1 R 1 

$G ORBITAL'/' ' ,T2R, ,T37, '-L- TIMp FLUX P.STML SnpjPlB^ 

SFLUX FLUX(SUM) •/' (OEG) (OPr,) (KM) (GAUSS) (P.R.) ( MRS ) Sl'F J ul H 3 

S #/CM*X=2/SEC ' ) POP I PI R4 

WRITE(6,202 )PSNTM1,PSMLN1,PSM.T1 ,PS.MAL1 »P.SNB) , PSMLl SnFTPlPP 

202 FORMAT!* ' » T^^l , F9 . 5 , T 2 , F7 . 2 . 1 X , F6 . ? , 1 X , FB . 1 , 1 X , F « . P , 1 X , F B . ? , T B(! , SOFTpI’^B 


S7!2X, 1PE10.3) ) 

210 IF!MOD! IPRINT,KPRIMT) .NE.O) GO TO 2?n 

TIFLUX = FLUXES !MRGYL\/)=-FLOAT!KPSTEP) -60. 


'>pp T Pi « 7 
p n F I p 1 p P. 
S 0 F 1 P 1 M P 


W R I T E ! 6 , 2 0 2 ) P S N T I M , P S N L 0 P S N L A T , P S M A L T , P S N B , P S M L Y 
SFLUX ES ! MRGYLV ) , T I FLUX , F L X SUf'' 

220 IPRJi\iT = IPRlNT + l 

3^ Sjc a{o{C ^ sjc ^ 0 K H I T L Z 0 F 1^ R ir" ^ K 0 fl 1'' ^ M f”) n 1 I }_ p 'i' 'I' n' ?|' 3|C 3;c 

THIS 'MGDULF N * U S T he used with RFRCFMT TF'^F h'' pfiiiLF 
STORE TIME IN INNER E OUTER ZONE^ F XTER I. 

IF ! PSML.lt .0.0. OR. PSML.GT. 11. '0) G(i TO ?Bn 
IZ = I ZONE! IMT! PSML/. 1 ) ) 


Snp T Pi 9i) 
SOFT Pi 9l 
snp I Pi 92 
-Is ;;s ;)£;[; ijc S I ’ F I P 1 9 3 

TD194 

i 0 P T P 1 9 3 
SHF ! Pi 96 
A n F ? p 1 9 7 


LCOUMTIIZ) = LCnUMTIIZ) + 1 
GO TO 260 

250 LCGUi\iT!4) = LC0UMT!4) + 1 
260 CONTINUE 

EXPOSURE INOEX MODULE 

STORE FLUXES AND TIMES F'i INTENSITY RANGES 
GO TO! 270, 2B0) , I PASS 
270 ISWTCH=ISWTCH+1 


3 0 F T P 1 9P 
sop T PI 99 
SilF T P200 
SUFI 9201 
;;c :;c ;;c :jc ;[c pi F | P 2 0 2 

SQF T P20 3 
SOF 1 P204 
SOF I P20 5 


280 INTRNG = ! 8- IN T ! 1 . 0-S I GN ( 0 . 5 , ALG FLX ( MRGYL V ) -7 . 0 ) ) 
$!7-INT!ALGFLX!NRGYLV) ) ) ) I NT ! 1 . 0 + S I GM ! 0. 5 , 
SFLUXESINRGYLV) -1.0009) )+l 


SOF I P2()6 
SOF I P207 
SOF T P2()« 


EXPFLX! INTRNG ) = EXPFLX! I N TRNG ) +F LUX E S ! NR GYL V ) -60 . F I. OA T ! SnFIP20 9 

SKPSTEP) S0FTP21O 


EXPFLX ! 10 ) =EX PFLX ! 10 )+FLUXES ! NRGYLV ) -60 . 0-FLOAT ! KPSTFP ) 
EXPTIM! INTRNG) = 6XPT I M ( INTRNG ) + FLOAT ! KPSTFP ) - .0166667 
EXPTIMdO) = EXPTIMIIO) + FLOAT ! KPSTFP ) =i- .0) 66667 

PEAK AND TOTALS PER ORBIT MOOULF 
DETERMINE ORBIT NUMBER AND TOTAL FLUXES PER ORBIT 
IFIPSNTIM.LT.TAU) GO TO 300 
PEAK = -1.0 

TAUFLX ! M0R8IT ) = FL X SUM-O F L X SM 

□ELXSM = FLXSUM 

NRBITO=NORBIT 

N0R8IT = MORBIT+1 

TAU = NORBIT PERIOD 

IF (NOR BIT. LE. 50 ) GO TO 300 


S0FIP211 
SOF I P2 12 
SOF I P21 3 
5^ n' 0 n F I p 2 14 
S0FJP21 5 
SOFl P216 
SOF I P21 7 
SOF I P2 1 B 
SOF I P2 19 
SOF I P220 
SOF T P221 
SOF I P222 
SnFIP223 
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301 

't' 

300 


41(' 
44 0 
4 S U 

0 


4 6 u 


47 n 
401 


NOkRIT EXCEcOS limit of 50, I ) 


ORBIT 


310 
C 'i' 

0 5 Ji 5 {? 5 jc 

0 'i^ '*' 


0 

0 'I* 

6.00 


WRITE(0f301) SDFIP226- 

FORMAT (■ OERkOR: MORBIT EXCEEDS LIMIT OF 50. > ) S0FIP225 

STOP SQFIP226 

DFTFRMIOE FLUX PEAKS AMD POSITIOMS PER ORBIT S0FTP227 

IF(FI.UXES(NR(;YLV).LE.PEAK) go to 31^ S0FIP22B 

PKFLX (MflRBIT ) = FLUXES ( NRGYLV ) SnFIP229 

RKTIMIMORBIT ) = PSNTIM S0FIP230 

PKLIIM ( MORHIT ) = PSNLON S0FIP2 31 

PKLATIMOKRIT ) = PSNLAT S0FIP232 

PKALTIMORBIT) = PSNALT .S0FIP233 

PKB(MORBIT) = PSMB S0FIP234 

PKL(R'ORBIT) = PSfvL S0FTP235 

PPA< = FLUXES! OR'^YLV ) S0FIP236 

CnOTIMUE S0FIP237 

^ ^ *|C 9|C S{0|C 3|C !|C ?!« f T r .* Nl A G M E T I C S B I P f I M ! ) D 0 I V V ^ ^1* 'i^ ¥ ^1^ 5,4 5j4 SjS 3jC ),C 3{C ^ I P ^ 2J P 

>;4 5;=S|4 4|4S,4;;=,4 THIS MODULE MUST RE USFD WITH SOLAR PROTON MODULE 5l=>i=>!'>;'5|«5|*S0F I P239 
I F ( INT ( PSNL ) .GR . S „nR .PSNL . LE . o . r«) L = L+ 1 SOFT P240 

5,4 5^ sl4 5l4 3[4 3[4 >|4 5jc sjs Sj4 5;c BLOCK 4 • L 0 U P I 0 G ( R t? A D“ L 0 0 P ENDS HERE) 5 J 4 '[4 3J4 3|4 3|4 5 I 4 3|4 3[4 3j4 3J4 3J4 ^ Q p J P 2 4 X 

IPASS=2 S0FIP242 

Gil TO 50 S0FIP243 

5i4 514 3 J 4 3j4 3l4 Sj4 3 J 4 3|4 3j4 5[4 3|4 3;4 3|4 3j4 ;J4 3 J 4 3 J 4 3p B L C 1 ( / K 5 * 0 U T P 1 j T P R E P A R A T 1 0 M 5^ 3i? 3;4 3j4 3J4 3J4 3j4 5J4 3;4 3^4 3^4 3j4 3j4 5(1 5^ 3^ 3)4 S Q E I P 2 4 4 

COMPOSITE ORBIT SPECTRUM S0FIP245 

AFCTRS = (KPSTeP*1440.0) / 

00 410 NRG =1,30 

AIFLXS(NRG) = AI FLXS ( NRG ) 5KAFCTR S 


( PSMT I f.-.<864no. 0 ) 


' 'r 'I' 't' 'I' 'I' 'I 


IF(AIFLXS!NRG).LE.0,0) go TO 440 
ALHFLX(MRG) = AL0G( AIFLXS(NRG) ) 

COMTIMUF 
no 450 MRG=l,29 

niFFLX(NRG) = AJFLXS(MRG)-AIFLXS(NRG+1 ) 
niEFLX(30) = AIFLXS(30) 

V 3,4 3,4 3 J 4 3 J 4 3 p 3{4 3,4 3,4 3,4 3 ^ 3|4 3^4 3,4 3,4 5 j' 3|4 3 [t 3 j 4 3|4 3)4 3 J 4 P ITT R C E N T T I H M 0 i ) t 1 L P 

THIS MOnULfc MUST BF USED WITH ORBIT L-ZONE BREAKDOWN MODULE 
CALCULATE AMD PRIMT PERCENT TIME TABLE 
LSUM = lCOUNT ( 1 ) + LCOUNT(.?) + LCGUNT(3) + LC0IINT(4) 

IFILSUB'.EQ.O) GO TO 4 70 
on 4H0 IL=1,4 

PTIME( I L ) =FLOAT ( LCOONT ( I L ) =i=KPSTEP ) . B6B67/TML AST . 

PT IZ = PTIME ( 1 ) + >JT IMF (?) 

WHITE! 0,401 ) PTI 7 , ( PT If'E ( I I ) , I 1 = 1 , 4 ) ,PSNTIr’ 

COwTINOE 

FORMAT (' 05.'":45X5;:5;4 ppkqrNT C>F TOTAL LIFETIME SPENT IMSIDB 
strapped particle RAOIATION BELT ^xcxcsici //6X , • INNER ^ 70NF 
■S2.8) ; ',F6.2,' /X8X, 'OUTSIDE TRAPPING RFOTOP' (1.0 

s: X '/IBX, ' INSIOE TRAPPING RROinM (T.l <= L < 


AND nijTSlDP 
( 1.0 <= L < 
<= L < 1.1 ) 
? . « ) : . ' , • 


SOFI P246 
S0FIP247 
S0FIP248 

SOF I P249 
S0FTP250 
S0FIP251 
SOFT P252 
SOF IP253 
5npTP254 
SnFIP255 
S0FIP256 
■SOF I P237 
SOFI P25« 
■SnFIPZSP 
SOF I P260 
SOF IP261 
SOFI P262 
SOF IP2B3 
SOFI P2/S4 
SOF I P2/SP 
SOF I P266 
S0FIP2B7 
SOF 1 P26R 
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' /6X, 'UUTEK ZONF (?.B <= I. <= 11 . 0 ) : s; ' / ft X , ' F X TESnF I P269 

■SRivAi. (L > 11.0) : '.Fft.?,' ^'//' TOTAL ORBIT TIME T b : • , SOF I P2 70 

‘^Fm.;?,' HOURS') S0FIP271 

Q p;i 1 F F E R E N T I A L SPECTRUM i-IOnULF =? >? 5;= -t- =1'. ^ >!: « « Jjt S 0 F I P 2 7 2 

CALL nSPCTR ( ALiMFLX( 1 ), ENERGY! I , ITYPE ), ni FSPC( 1 ) ) _SnFIP273 

'• j ¥ V V 5|C 3|C S[! S|5 ?|C 5iJ 5|^ '|C 3|C Sjf ^|t 5,C 2,% S t H__ A R PROTON 0 0 U ^ F ^ 'I> ^ 'i^ ^ ^ ^ 'i? 'JS Sjc ^ 3[5 >|S ^ 0 F I P 2 7 A' 


C This MODiJLE MUST BE USFO WITH GEOMAG. SHIELOING MonULF 

T=l?. 

IT = T 

10=90 

ISWTCH=IS'.a)TCH + 7 
IF(L.LF.O) go to 510 
CALL SOLPROIT , 10, F, INALS ) 

EXPUTM=FLOAT ( L=i=K P ST E P ) =:= , 0 1ft ftft ft7 


^=«*5;=SnFIP275 
SOFI P276 
S0FIP277 
SOFI P2 7ft 
S0FIP279 
SOFI P280 
S0FIP2ftl 
S0FIP2ft2 


C 

C 


0 ! ) 

^ 1 ' ^ 


'r* ^ 'f* 


EXPFCT= ( EXPOTM/PSNT IN ) 
no 500 J = l,2i./ 

F ( J ) =F { J ) =:<EXPFCT 


SOF I P2ft3 
SOFI P2ft4 
S0FIP285 


CONTINUE S0FIP2R6 

CONTINUE S0FIP2B7 

OUTPUT PUNCH MOOULF =1- s;-- sjt :|= W >!i >1: 5|< 5|! S 0 F I P2ftft 

POOCHES ENERGY, IMTEG AMO OIFF FLUX, SOLAR PROTONS IF PRESFMT S0FIP289 

WRITE! 7f 605) n anE , I MCI , I PRO . I APG , MODEL , BLT I ME S0FIP290 

WRITE ( 7,ft02 ) ( ( ENERGY! ( I I-l ) 'i'ft + JJ, I TYPE ) , JJ=1 , ft ) , XI.AftFL ( 1 , ITYPF ) , I ISOF IP291 


■^,11 = 1,5) 


S0FIP292 


WR ITE ! 7,602 ) ! ! AI fLXS ! ! I I-l ) -6 + JJ ) , J.)=l ,6).XLAREL!7,ITYPF),II,II = 1,S0FIP2 93 
s5 ) S0FIP294 

WR I TE ! 7 ,602 ) ! ! DI FSPC ! ! I I - 1.) *ft + J J ),JJ = },ft),XLA8FL!3,ITYPF),II,II = l,S0FIP295 
65) ^ S0FIP296 

IF!L.LF.0) go to 600 S0FIP297 

WRITE! 7,603) Ii.),T, O' ALF,EXPFCT S0FIP298 

WRITE ! 7, ft04 ) ! ! SPNRG! ! I I-l )^!=5 + JJ ) , .1 J=l, 5 ) , PR0TL8 ! 1 ) , I I, I 1 = 1,4 ) S0FIP299 

WRITE! 7,604 ) ! ! F ! ! I I - 1 ) - 5 + J J ) , J J= 1 , 5 ) , PR OT LB ! ? ) , I I , 1 1 = 1 , A ) SOF I P300 


600 comtinue 

602 FORMAT! 1P6E12.4, Aft, I?) 

603 FORMAT !' SOLAR PROTONS aENERGIES=?0 0=',I2,' TAU=',F4.1, 


S0FIP301 

S0FIP302 

S0FIP303 


ft' N ale =',!!,' EXPFCTR= ' , F5. 2 ) S0FIP304 

ft'Jft- F0RMAT!1P5E12.4,I2X,A6, 12) S0FIP3O5 

605 format ! 3A4, IX, 12 , ' / ' , I 5 , , 16, IX, ' I ! a/CM=i=':=2-SFC ) D! /CM=.-*2-SEC-KES0F I P306 


ft V ) MOO/TM= ',11, ' / ' , Fft . 1 ) 

C If; sic OUTPUT TABLES MODULE 1 

DO 900 NTBL=1 ,NTABLS 

WRITE!ft,2) MODEL, MODLRL,BLTJMF, NAME, INCL, IPRG, 
SXAMwIN 

GO Tt i ! 7 1 0 , 7 ' >0 , 7 3 O’, 7 2 :' ) , I ,S W T C H 
C '*** COMPOSITE ORBIT SPECTRUM ANO EXPOSURE INDEX 


S0FIP307 

s|c s|c sjc >|c s|c s|c s|c s|c s|c s|c s|c s|c sjc s|c s|c s|c s|c s|c s|c sjc S 0 F I P 3 0 8 

S0FIP309 

I APG, ITAPE,PPRIOD, S0FIP310 

S0FIP311 

S0FIP312 

S0FIP313 
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700 WRITE (6,701) ( T YPL BL ( K , I T YP E ) , K= 1 , 3 ) , EMEROY ( MROYI. V , I T YPF ) . 6nFIP314 

S (ENERGY (N, I TYPE ) , AI Ft. XS ( N ) , f) I FFLX ( W ) , n I F S PC ( N ) , F I RNGS ( ^' ) , 6PFT P31 6 

SFIRNGS (N + 1 ) ,EXPTIM(N) ,EXPFLX ( M ) , M=1 , 10) , ( FVERr,Y( P, TTYPF) . A ? FLX\S ( M ) ^riF T P316 
S , D I FFLX ( N ) , DI FSPC ( N ) , N= 1 1 , 30 ) .Snp I P31 7 

701 FORMAT ( ' + • , A 1 X , 16 ( ' # ' ) , 3X , 3 A4, ?X , 16 ( ' ' ) / ' 41 X , A9 ('=?•)////' li ' , SHPIPblH 

ST18,15( '*' ), ' COMPOSITE ORBIT SPECTRUM 1 6 (•>:'•), T«0 F X PPSi !R c ^0= I p31 9 

SIMDEX:ENERGY>' ,G9.2,Tll?. 'MEV '/' 0 T 1 B ENERGY AVFpAGEO 90FIP320 

s DIFFERENCE AVERAGED DIFFE- . IMTEN^ITY FXPOSURP TpS0FTP321 

STAL « OF*/' ' ,T18, 'LEVELS INTEGRAL FLUX l-MTEGRAL FLUX P FI'- I 1 SOF I P322 

3.AL FLUX ', lOX, 'RANGES DURATION ACCUMUL ATFO ' / ' ' , T1 «, ' >( ■'FV' ) .snFIP323 

S #/CM>J'*2/S6C *4/CM“'’i'2/SEC/DE a /CM-=i' ? / SEC /■< F V a/C’i*=.'2/SFCS0F I P326 

$ (HOURS) PARTICLES'/'O' ,T1B,0PG9.4,T?3, ' ',IPF9.3,7X, SnFTP325 

S1PE9.3,8X, 1PF9.3,T81 , 2AA,TR1. , * ZERD F LUX ' , 1 X , OP FI 0. 3 , 1 X , 1 PR 1 3 . 3/ « ( ' SOF I P326 
S ',T18,0PG9.4,1PE9.3,7X,1PE9.3,8X,1PF9.b,TB 1, A4,»-',A4,i.x,OPFl0.3,SnFTP3-2 7 
S1X,1PE13.3/),' ' ,T18,0P(G9.4, IPE9.3, 7X , 1P69.3, RX, 1 PF9.3 / ‘ ' ,T8l , 24 4S0F T P3/<« 


S,T81, ' 


TOTAL ' , IXtOPFlO. 3, 1 X , 1 PF 1 3 . 3 , ?P- ( T 1 8 . 0 PG9 . 4 , 1 PE9.3, 7X, 


I 


) ) 


ITYPE) ,/'IFLXS(i'' ) 


S1P69.3,8X,1PF9.3/' 

GO TO 750 

C =«'=«'* COMPOSITE ORBIT SPECTRUM ONLY 

710 WRITE (6,702) ( T Y PLBL ( K , I T YPF ) , K= I , 3 ) . ( FM FR GY ( i 
sniFFLX(M).,DIFSPC(N) ,N=] ,30) 

702 format ( *-P ' ,41Xf 16( '« ' ) , 3X. 3A4, ?X, ] 6 ( ) / • ', 4;: x, 49( '«=')////' O ' . 

ST40, 15( '-!='),' COMPOSITE ORkIT SPECTRUM ', 1 . 6 ('-')/' 0 ' , T40, ' fnfrGY 
S AVERAGED DIFFERENCE AVFRAGFO DIFFF-'/' ' , T40 , ' I. FVFLS 

s INTEGRAL FLUX INTEGRAL FLUX RFNTIAL FLUX'/' ' , T4() , • >( 'V c\/ ) 

$ #/CM**2/SEC a/CM>i'* 2 /SFC/nF a/C.M>i'='''?/SFC/KFV' //30( ' ',T4l. 

iOPG9.4,T46, ' ',1PE9.3,6X,1PF9.3,HX,1PE9.3/)) 

GO TO 750 

C COMPOSITE ORBIT SPECTRUM WITH SOLAR PROTONS AMD EXPOSURE INOEX 

720 WRITE(6,703) ( TYPLBL ( K, ITYPE ) ,K=1 ,3 ) , IT, 10, I N AL F , EMcr GY ( t'R GYL V , 

%I TYPE ) ,EXPFCT, ( ENERGY! N, ITYPE ) , A 1 FLX S ( N ) , 0 I FF L X ( M ) , D T FSPC, ( M ) , S PNPGS OF I P 344 
5K (N ) ,E ( N ) , FIRNGS ( N ) , F IRNGS ( N+1 ) , EXP T I M ( N ) , F X P FL X ( N ) , w= 1 , 1 0 ) , ( Fimppp- YSOF I p345 
$(N, ITYPE ) ,AI FLXS (N) ,01 FFLX ( M ) ,DIFSPC( N) , SPMRG( N) , F ( N ) , N= 1 1 , 20 ) , SUFI p346 
S (ENERGY (N, ITYPE ) , AIFLXS (N ) , 01 FFLX ( N ) , HI FSPC ( N ) , N=21 , 30 ) S0FTP347 

703 FORMAT! '-P' ,41X ,16( '=s') ,3X,3A4, 2X, 16( '-')/ ' ' , 41 X , 49 ( i ^ i ) / / / / i O i , ,sOFIP34H 

S62X, '***=:= SOLAR PROTONS >:==.'❖=!='/ /63X ,' FOR T AU= ' , I 2 , ' , 0= ' , I 2 , ' : NaLP = SGF I P349 
$' ,I 1/3X, 15( '*' ) , ' COMPOSITE ORBIT SPECTRUM 1 5 ('-'), 5X ,' m T TM GEnMS0FTP350 


SnFTP32.G 
Sf’F T P330 
,'>0FTP331 
^DF T P 332 

S''FIP333 
SOFT P3i4 
SOr )p335 
SOF I P3i6 
SnFTP337 
•SUE I P338 
SQF 1P339 
SUE I P 3 41.1 
SOF I P34I 
SUE I R342 
SQF IP343 


$AG SHIELDING' ,8X, EXPOSURE INDEX: FNER GY> ' , G9 . 4 , T 1 ? 5 , ' OFV 


S64X , ' ( FXPOSR FACTOR= ' , F4. 2, ' ) ' //3X, ' ENERGY 


SEREMCE 
$SITY 
»L FLUX 
$S 

sSbC/DE 

•fi/SEC 


AVERAGED DIFF- ENERGY 

EXPOSURE TOTAL * OF ' /3X LEVELS 
RENTIAL FLUX LEVELS 

DURATION ACCUMULATED' /3X, 'X MFV) 
«/CM=l--=l=2/SEC/KEV XNEV) 


ayFRAGrO 

total 

INTEGRAL PLUX 
FLUFNCE 

-a/CN=i=^2/SFC 

a/c\-Px>X2 


' /SOF IP3sl 
niFFSOFl P3d 2 
INlENSriF I P353 
TWTEGRAS0FIP354 
R AnpfSOF I m355 
i/CN=?Y2./S0FT D356 
.:/CM!!==;c2SnF I P357 


(HOURS) 


PARTICLES' //T4,0PG9. 4, T9, ' 


, 1PF9.3,6X, 


SOF T P3 5B 
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$1PE9. 3,7X, 1PE9, 3, 1 1X,0PF4. n,7X, 1.PF9. 3, T95, T9'5, ' 7-KO FUlX' , 

£0PF11.3,'lPE14,3/8( T4,0PG9.4, T9, ' ' , I P ?=9 . 3 , 4 X , I PP 9 . 3 , 7 X , 1 m H9 

$ , 11X,0PF4.0,7X, 1PE9.3,9X, A4, ' , A4, OPFU. . 3, 1 PE14.3 / ) ,T4,0Pf;9.A, -j q 

S' ' ,1PE9.3,6X, 1PE9\3,7X, 1PE9.3, 1 IX , OPF^. n, 7X , I PF9. 3 , T93 , ^o4 , 

ST95, ' TOTAL ' ,OPFl 1 . 3 , 1PE14. 3/10 ( T^, ''''PG9.4, T9, ' ' , 1 pp9. 3 , 4 

$,1PE9.3, 7X,1PF9.3, 11X,OPF4.0,7X, 1PE9.3/) , VM T4,0Pf;9.4,TP, ' ' 

S1PF9.3,6X,1PE9.3,7X, 1PE9,3/) ) 

GO TO 750 

C COMPOSITE ORBIT SPECTRUM WITH SOLAR PROTONS 

730 WRITF(6,704) (TYPLBLI k, ITYPE) ,K=1,3) , IT, lo, IMA'.?, EXPECT, ( FHPROV ( , 
SITYPE) .AIFLXS(M) ,DIFFLX(K') ,niFSPC(M) , SPORGI '-i ) ,F(m) ,0=} , ?0 ) , 

S( ENERGY ( M, I TYPE ) , AIFLXSI N ) , D I F F I. X ( N ) , n I FS PC ( f- ) , 0 = ? 1 , 3(0 
704 FORMAT ( ' + ' ,41X, 16( ) , 3X , 3A4, ?X , 14 ( ' - ' ) / ' ' , 41 x,4q( • =? ' ) ////T93 , 

SOLAR PROTONS xc>;=5,''-:k ' / /T93 , ' FOR T Al != • , I ? , ' , n= ' , 1 P , ' : W«LE=', 
SI1/19X,15( •#•),» COMPOSITE ORBIT SPECTRUM ', 1 5 {•»•), 1 HX , ' wTTw Gpn 
SAG SHIELDING '/T94, •( EXPOSR F ACTOR =', F4 . p ,')'// 1 9 X FNFR C-Y A\/E 

SAGED DIFFERENCE AVERAGED n T F F- ' , P:'.X , ' EN p R GY TOTAL 

S/ 19X ,' LEVELS integral FUIX INTEGRAL FLUX REOJIAL FLOX',/PX, 


S 'LEVELS 


FLOENCE ' /19X, 'XMEV) 


Jt / E f.'i ? / s Ff C. S / C ‘R * ? / S ^ C / I" E 
S #/CM##2/SEC/XEV» ,20X, 'XMEV) » /C Msecs' 2 • / / ? 0 ( T P 0 , 0 POg . 4 , TP S , ' 

S ',1PE9.3,6X,1PE9.3,7X, 1 PE 9 , 3 , P4X , ' 'P p4. o , 7 X , I P P 9 . 3 / ) , 1 0 ( TPO, 

S0PG9.4,T25 , ' ' , 1 P E9 . 3 , 4X , 1 P F 9 . 3 , 7 X , 1 P E9 . 3 / ) ) 

750 CONTINUE 


Q 5i< Jjc 3^ ^ 


' 5jC 5{S 5{i 5JC SjS 


'i* -V 'C- 'i" '1" ' 


OUTPUT TfiRLFS '^OnULP ? 


S'5F T P 
.•5S0F T 
, SHF I P3'Sl 
P36^ 
VSHF IPB'Sji 
, SOP T 
SnPTP3A5 
SOr T P366 
P3.S7 

SOP T 

S ^ ; F I P 3 ^ ^ 
SOFT P3 70 
SPFTP371 
SOF I P3 7/^ 

; SnFTP373 
PsilP] P374 
» SHF I P37b 
STm- T P37R 
SOF 1^^377 
SPF T P37« 
SOP I P37P 

Spe I p3p() 

SHF TP3^U 
=1' S n F T p 3 p 3 


I PRG> I Apr;, TTAPPt PF'RIOOi 


S0FT^3p3 
SOF I 


C PFAK AND total FLUXES PER PERIOn 

WRITE(6,2) HOD EL t MOO L HI. , BLT I ME , NAME» I NCL 
5XAMM1M 

WRITE(6,801)(TYPLRL(K,ITYPE),K=1,3), ENERGY! WROYLV, ITYPF ) , SOFI P3«4 

S (N , ( PKVALU(M,K ) ,K=1 , 8 ) , M=1 , NRPI TO ) SnFTp3»7 

801 FORMAT! ' + ' , T35 , ?4 ! ' * ' ) , 3X , 3 A4 , ? X , ? 7 ! ' -^ ' ) / ' ',T35, SnFlP3«B 

TABLE OF PEAK AND TOTAL FLiJXES PER PERIOD : ENERGY > ' , <;o . 2 , T9760F I P 3« 9 

S, 'MEV *>s'/' ' , T35 ,48 ! ' ' ) // ' 0 ' , 13X , ' PER tod PEAK FLUX POS0FTP390 

SSITION AT WHICH ENCOUNTERED ORBIT TIME FIELD! B) LIN'=!L) SDFIP3P1 

S TOTAL FLUX'/' ', 13X ,' NUMBER EMCOUNTERFD LONGITUDE LAT T TODSOF I p392 
SE ALTITUDE', 41X, 'PER ORBIT'/ ' ' , ?3X , ' #/CM-”-?/SFC ' , ? ! 5X , ' ! DF,- ) ' SOF I P3Q3 

S) ,6X, ' !KM) ' ,7X, ' !HOURS) ' ,6X, ' !GAUSS) ! F.R. ) & /CM=X=? ?/nR« 1 I 'S0FTP394 

S// ! ' ' , 14X, 14, 1PE14. 3,0PF13. 3, FIO. 2, FI ?. ?, F13. 5, FI 7 , 5, FIO. 7, 1PF15.S0F 1P395 

S3)) S0FTP394 

C *#!;! X! xs xs Sr BLOCK 6: PROGRAMi TERM I NAT I ON xsxsx^x^xtx^xtXtxtXi^fSrJrXorXsxiSOF I p3_97 


900 CONTINUE 
GO TO 10 
999 STOP 
END 


SOF I P3QB 
SOF I P399 
SOFI P400 
SOF I P401 



Q ^ ^ ijc 

f , ^ 

^ 'I' 

C, 5|< 5jc 5,< 

Q 5;< 5|i 5|C 

(^ jJC i\i 5? 

Q ?[c 5jc ?JC 

(, 

Q >? s;t 3jc 

f , 5;i 5jc 


C*, 'i' 'I' 'I' 


C 'J« 

Q 5j; jJh 


1='^ 

15 

Q 5|< s’? 

(3 'I' '1^ 


30 


40 


sjc s^ s;c S.C ^ s|? s;c si? s|c sjc s{c s|? s!? D I F F F 1^ h" N T I ^ L S P F 0 T K 1 1 M S FJ HKOUTI Ntr ^^''^'<^^^^^^^5{t3jC5jcsJc5jc^Qp^pQ0O2 
CALCULATES FIRST OERTVATIVES OF INPUT SPECTRUM OFFINFO BY FF VS XXDFSPC004 
I'mPUTs XX ~ 3(> INTECRiiL THRESHOLD ENFRCIFSf IN MFV ( R=i^4 ) DFS PC006 

FF - ALOG OF THE INTEGRAL FLUXES FOR THE. 30 ENERGY ( R>J^4 ) DPS PC008 
LEVELS, IN PARTICLES/CN^^'?/SFC DFSPCOlO 

OUTPUT: OD - differential FLUXES OBTAINED FRON TMF INTEGRAL ( R>!^4 ) DFSPCO 12 

FLUXES, liM FARTICLES/Cf-'»?>i'?/,S5C/KEV DFSPC014 

^ 5,C ^ ?}c 5JC 5}c 5*|C 5|C 5jc ^ 5JC 5|C I»JC ^ 5jC 5jC ?jc ?[C jjC 5^ 5|€ ^ 3|C ^JC ^|C ^ ?JC ^ ^ 5|% ^jC ^ ^ JfjC ?|C ^ ^ ^JC Jjt 3^ 3{C ^ ^ S P 0 X 

THIS IS A MODIFIED VERSIOr'i OF A PROGRAM (DfSIFIJ) OBTAIMED FROM DFSPCOlft 
IMSL LIBRARY 1: fiUT HOK / I MPLEM ENTOR - G.L. SMITH DFSPC020 


:(C 5|s # i;c i|« >'.5 # Jis # 5,t >;£:<« si! * !|s ;Jt j;s 5? Jit j;« « SU 3S **** =l« 5>; »! >r V >;«>!!>?: =i5 5? !i! 5|S Sli =1< # Jit 5|« ^ V 5|S j|< 5i: !f: ❖ 5i= sis # s!« D F S P C 0 2 2 


SUHRflUTIiME DSPCTkI FF.XX.DD) 

I'-'PLICIT KEAL=«^H ( A-M.O-? ) 

REAL=^4 DD.FF.XX 

DIMRmSKjN F(3D) ,X(30) ,D(B0) ,H(500) ,FF(30) ,XX(30) ,DD(30) 

DATA EPSLN , GM6GA/1 . n-ft , i . n? 1796 HDH/ 

DATA INITIALIZATION 
M = U 

DO 3 L=l,30 
nn(L)=o.o 

DcTERMINE SIZE OF ARRAY: OBTAIN M S K INDICES 

M = ri OF NONZERO FLUXES - l: K = » OF NONZERO FLUXES 

DU 10 K=l,30 

IF(FF(K) .eo.O. ) GO TO 15 

h=K-l 


DFSPC024 

0FSPC026 

DFSPC028 

DFSPC030 

DFSPC032 

0FSPC034 

DFSPC036 

DFSPC038 

DFSPC040 

DFSPC042 

DFSPC044 

DFSPC046 

DFSPC048 

DFSPC050 


FK)=FF(X)+AL0G( 1000. ) 

X(K)=XX (K )=itioOO.DO 
D(0=X(i<) 

K=M+1 

IMK.LT.IO) GO TO 170 
SMOOTHING INTEGRAL FLUX 
CALL SMOOTH! X, F.M ) 

CALCULATE SECOND DERIVATIVES USING CENTRAL niFFFRFNCES 
DU 30 I=1,M 

H( I )=X ( I+l )-X( I ) 

H(K +I ) = ( F( I + l )-F( I ) ) /H< T) 

DO 40 

H(?.Ji=X + I )=H( I-l )+H( I ) 

H(3*K+I )=.5*H( I-n/H (?*K+I) 

H(4=(=K + I ) = (H(K+I )-H(K+I-l) )/H( ?=itK+I ) 

H(5*K+I )=H(4>:tK + I )+H(4>i^K+I ) 

H(6*K+I )=H(5*K+I )+H(4JitK+I ) 


0FSPC052 
0FSPC054 
DFSPC056 
OFSPCasa 
DFSPC060 
0FSPC0G2 
NFS PC 044 
DFSPC066 
DFSPC068 
DFSPC070 
DFSPC072 
0FSPC074 
DFSPC076 
DFSPC078 
DFSPC080 
DFSPC082 
DFSPC0R4 


H(5=!=K + 1 )=0. 


0FSPC084 
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H ( 6»I=K ) =0. 


DFSPC0R8 


^ Jjt S|5 

*5 0 


60 

70 

Q ^ ais ;,x 
80 


9') 
1 ()0 


BFGIM ITERATION ON SECONO DERIVATIVES 

Kr.nuNT=o 

ETA=0. 

KCnUOT=<CmiNT+l 
DO ,7 1) I=?,M 

w= ( H ( 6>!‘K+ I ) -H( 3^'K+I) *H ( 8*K+I-1 )-( . 5- 
S I) IOMEGA) 

IF (DABS(W) .LE.ETA) GO TO 60 . 
FTA=DABS(W) 

H ('b ’I«|< + 1 ) = H ( 5 =:= K + I ) +W 
COOT I MO E 

I F ( KCOOMT.GT .5*K )G0 TO 170 
IE (ETA.GE.EPSLM) DO TO 50 
CniVVERGEMCE OHTAIMED 
DO 80 I=1,M 

H(7’:'K + I ) = (H(B=('K+1+I )-H(5*K + I) )/H( 1 ) 

Dl l I AO J = 1,K 
1 = 1 

IF (D( J).ED.X( 1) )G0 TO 130 
IF (D(J)-X(K )) l'i0,110,110 

IF (D(J)-X(I)) l?n,l3^,100 

1 = 1+1 
GO TO 90 


0FSPC090 

DFSPC092 

DFSPC09A 

DFSPC096 

DFSPC09B 

H( 3*K+I ) ) AH( 5=SK+ T + 1 )-H( 5#K + 0FSPC100 

DFSpCIO'2 

OFSPCIOA 

DFSPC106 

DFSPCIOS 

DFSPCllO 

DFSPC112 

DFSPCllA 

DFSPC116 

DFSPC118 

DFSPC120 

0FSPC122 

DFSPC12A 

0FSPC126 

0FSPC12B 

DFSPC130 

UFSPC132 

DFSPC13A 


110 I=K DFSPC136 

12u 1=1-1 DFSPC138 


C --- COMPUTE D( J) DFSPC140 

130 HTl = n ( J >-X ( I ) DFSPC1A2 

HT? = D ( J )-X( 1 + 1 ) DFSPCIAA 

PPnn=HTlAHTP 0FSPC1A6 

H ( 8AK+J ) =H( 5AK+1) +HT1*0 ( 7*K+I ) OFSPClAp 

nFLS0S = (H(5*K + I )+H(5AK+l + I )+H(8=!^K + J) )/6. DFSpC150 

1^') D(J)=-(H(K +I )+(HTi+HT?)ADELSOS+PR0DAH( 7+K+I )-. 1666667) DFSPC152 

C SNOOTHIMG DIFFERENTIAL FLUX DfSPC154 

call SMOOTH! X,D,M) 0FSPC156 

on 160 1=1, K DFSPC158 

F ( I )=? .7182818?8D0**( F ( I )-ALOG( 1000. ) ) DFSPC160 

160 DD(I) =D(I)^F(I) DFSPC162 

170 RETURN DFSPC16A 

EiviO DFSPC166 

DFSPC168 

+ SMOOTH DATA BY 3-POIMT AVERAGING OVER EOUAL INTERVALS DFSPC170 

SUBROUTINE SMOOTH ( X . F . M ) . 0FSPC172 

IMPLICIT REAL*8 ( A-H,0-Z ) 0FSPC17A 

OIMEIMSION X(30),F(30) 0FSPC176 
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i FINTER( XI ,X2, X3, Y 1, Y2, Y3, XIM )=Yl=i=( XTM-X? )=:M XIN-X3 ) / 

$( (X1-X2 )>i=(Xl-X3 ) ) + Y2=M XIN-Xl XI N-X3 ) /( ( X2-X1 X2-X3 ) ) 

$ + Y3-(XIfv-Xl )*(XLN|-X2) /( ( X3-X1 )=!=( X3-X?) ) 

= C 

i FI = F ( 1 ) 

■ DO 20 I=2,M 

sizei = X ( I ) - X( I-l ) 

SIZ62 = X( I+l) - X( I) 

C CHECK FOR equal STEPSIZES 

IF(DABS(SIZE1-SIZE2) .LT. 0.001) 00 TO 200 
IF(SJZE2. GT. SIZEI ) GO TO 210 

C STEPSIZE DECREASES - FIT CURVE AMD IMTFRPOLATF RACKWARO 

F2 = F( I+l ) 

X inter = X( I ) - SIZE? 

FI = FINTER ( X( I-l ) , X( I ) ,X ( I+l ) , F I , F( I) f F?, XINTER ) 

GO TO 300 

C STEPSIZE INCREASES - FIT CURVE AND INTERPGI ATF FORWARD 

210 FI = FI 

X INTER = X( I) + SIZEI 

F2 = FINTERI X ( I-l ) , X ( I) , X ( I + l ) , FI , F( I ) , F( I + l ) , X INTER ) 

GO TO 300 

C STEPSIZES ARE EQUAL - AVERAGE OVER EXISTING VALUES 

200 FI = FI 

F2 = F ( I+l ) 

*=1'* PERFORM AVERAGING 
300 FNEW = (F1 + 2.0=^F( I )+F?) Z4. 

FI = F ( I ) 

F(I) = FNEW 
20 CONTINUE 
RETURN 
END 


DFSPC17R 
UF-St^LlFO 
DFSPC 1 R2 
OFSPCl M4 
OFSPClftO 
UFSPLl HR 
0FSPC190 
OFSPCl'^2 
0FSPC194 
L)FSPC19F 
OFSoClQR 
L'FSPC2U0 
0FSMC202 
UFSPC204 
DFSPC20(S 
UFSi^C2 0R 
0FSPC21 0 
OFSRC212 
0FSPC21 4 
l)FSPC216 
i)FSwC21 8 
UFSPC220 
UFSOC222 
0FSPC224 
0FSPC224 
OFSPC22R 
NFS pc 2 30 
0FSOC232 
0FSPC234 
DFSPC 2 36 
0FSPC2 3 8 
0FSPC240 
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SUBROUTINE SOLPRO(TAl), 10, F, INALE) SOLPROlO 

C *** MODIFIED 9/77 TO RETURN INALE(# OF AL EVENTS) TO CALLING PROGRAM S0LPR020 
C interplanetary solar proton flux at I AU (FROM e>10 TO E>200 MEV S0LPR030 

C *** FOR ANOMALOUSLY LARGE (AL) EVENTS AND FROM E>10 TO E>100 MEV FOR S0LPR040 
C *** ordinary (OR) EVENTS) S0LPR050 

C *** SINGLE PRECISION DECK IN STANDARD FORTRAN IV FOR IBM 3B0 MACHINES S0LPR060 
C *** (EBCDIC, 029 PUNCH) OR OTHER COMPATIBLE SYSTEMS, S0LPR070 

C program designed and tested by E.G. STASSINOPOULOS, code BOl, S0LPR080 
C NASA GODDARD SPACE FLIGHT CENTER, GREENB'ELT, MARYLAND 20771 . S0LPR090 

C jitjjc* jjoje ### its* =(5 Jis SOL PR 1 00 

c INPUT; TAU MISSION DURATION IN MONTHS (RFAL*4) SOLRRllO 

C ***# 10 CONFIDENCE LEVEL THAT CALCULATED FLUENCE F(N) S0LPR120 

C ***♦ WILL NOT RE EXCEEDED (INTEGER*^) S0LPR130 

C **** OUTPUT; FIN) SPECTRUM OF INTEGRAL SOLAR PROTON FLUENCE FOR S0LPR140 

C ENERGIES E>10*N (1=<N=<20) FOR AL EVENTS S0LPR150 

C ENERGIES E>10*N (1=<N=<10) FOR OR EVENTS S0LPR160 

C ***# INALE « OF AL EVENTS FOR GIVEN TAU AND 0 S0LPR170 

REAL NALEtNALECFI 7,20)/-.157l, .2707,-.1269E-l,.442BE-3,-.fllB5E-5, S0LPR180 

«.77S4E-7,-.2939E-9,-. 1870, . 1951,-. 6559E-2, .1990E-3,-. 341 8E-5, SOLPR190 

t. 3740E-7,-.1599E-9, -.2007, .1497, -.31 79E-2, ,5730E-4,-.4664E-6, S0LPR200 

$,1764E-8,0. ,-. 1882, .1228,-. 1936E-2, .26 AOE-4,-. 1022E-6, 2*0, , S0LPR210 

*-.2214, .1149, -.1871 E-2, .2695E-4,-. 11 16E-6, 2*0. ,-.2470,. 1062, SOL PR 220 

*-.1658E-2,.2367E-4,-.9465E-7,2*0. ,-.2509,.87l0E-l,-.B300E-3, S0LPR230 

S.8438E-5,3*0.,-.2923, .8932E-l,-.1023E-2, .1029E-4,3*0. ,-.3222, S0LPR240 

*.8648E-l,-.9992E-3, .9935E-5, 3*0. ,-.3518, ,841 7E-l,-.1000F-2, S0LPR250 

*. 9956E-5, 3*0., -.3698, .795 lE-1,-. 8983E-3, .8940E-5, 3*0. ,-.2771 , S0LPR260 

$.5473E-1,-. 1543E-4,4*0. , -. 28 1 8 , . 5072E- 1 , . 25 1 1 E-4, 4*0. ,-.2845, S0LPR270 

S.4717E-1, .5664E-4,4*0.,-. 2947, .4405E-1 , .8507E-4, 4*0.,-. 2923, S0LPR280 

$.4111E-1, .1106E-3,4*0. ,-.2981 ,.3853E-1 , . 1 3 1 2E-3 , 4*0. , -. 3002 , S0LPR290 

*.3585E-l,.1529E-3,4*0.,-.3001,.3312E-l,.1781E-3,4*0.,-.3141, S0LPR300 

*.3248E-1, ,1654E-3,4*0./,F(20) ,G(20 ) S0LPR310 

REAL ORFLXCI 5,9) /.154047E3,-.522258E4, .714275E5,-.432747E6, .955315S0LPR320 
*E6, . 198004E3,-.448788E4, .438 148 E5,-.196046E6, .3255 2E6, .5291 20F 3, S0LPR330 

*-. 122227E5,. 112869E6.-.465084E6, ;710572E6, .121141E4,-.266412E5, S0LPR340 

*,226778E6,-.85728E6, . 120444E7, .452062E4,-.103248E6,,896085E6, S0LPR350 

*-.346028E7, .499852E7, .272028E4,-.499088E5, .35305E6,-. 111929E7, S0LPR360 

*.133386E7,.275597E4,-.469718E5, .314729E6,-.960383E6, .1 1165F7, S0LPR370 

*.570997E4,-.799689E5 , . 38 1074E6, - . 610714E6 , 0. , . 101E3,4*0./ S0LPR380 

INTEGER INDEX(20) /2*7,6,3*5,5*4,9*3/ S0LPR390 

1 FORMAT!' TAU=',F4.0,' 10= ', 1 3 , 3X ,' PARAMETER ( S) EXCEED PROGRAM LIMIS0LPR400 

*TS') S0LPR410 

2 F0RMAT(2X, 'FOR THE COMBINATION OF TAU AND 10 GIVEN, NO S IGN I F IC ANTS0LPR420 
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% solar proton fluxes are to re expected* TAU=^SF6*?,* 10=*, 12) 

IF(TAU.GT*72**0R, I0.LT.80)GD TO 500 

IP=100-I0 

M=IMDEX( IP) 

nalf=o*. 

DO 300 J=1,M 

300 MALE = MALE+NALECF( J, IP)^^'TAU=!'^( J-1) 

I MALE=NALE+1 .0001 
IF(INALE.GT.O) GO TO 400 

c calculations for or-fvemt conditions 

I T=TAU 

TF( IT.EQ.l.AND.IP.GTel6) GO TO 700 
P=FL0AT( IP)/100. 

0F = 0. 

DO 100 J=l,5 

100 OF = OF + nRFLXC( J,IT)>!^ P^J^^ ( J- 1 ) 1 . E 7 
E = 10. 

DO 200 N=l,10 
G(N)=EXP( ,015R>J'(30.-E) ) 

F(N)=0F^G(N) 

200 F-E+TO, 

GO TO ROD 

c calculations for al-event conditions 

400 E=10, 

DO 600 N=l,20 

F(N)=7.9E9^EXP( ( 30.--E ) /26. 5) FINALE 
600 E=E+10, 

GO TO ROD 

700 WRITE(6,2) TAU, 10 
. GO TO ROO 

500 WRITE (6,1) TAU,IO 
ROD RETURN 
END 


S0LPR430 
S0LPR440 
S0LPR450 
S0LPR460 
S0LPR470 
S0LPR480 
S0LPR490 
SnLPR500 
S0LPR510 
S0LPR520 
S0LPR530 
SOL PR 540 
S0LPR550 
S0LPR560 
S0LPR570 
S0LPR5R0 
S0LPR590 
SOLPR600 
S0LPR610 
S0LPR620 
S0LPR630 
S0LPR640 
S0LPR650 
S0LPR660 
S0LPR670 
S0LPR680 
S0LPR690 
S0LPR700 
S0LPR710 
S0LPR720 
S0LPR730 
S0LPR740 
S0LPR750 
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■ 4 * SOFIP TEST ** 

♦* 90DBG/ 2000KM/ 2000KM ♦* 

4c ««««* 41 # 4c «* 4 ( 4 ^ 4(« 4c « 


T-1 

Header Page 





NAME = TEST 90/2000 
I NCL « 90 

I PRG - 2000 

lAPG = 2 0 00 

I TAPE « 5116 
MODEL = 5 

PERIOO= 2« 119969 
8LTIME* 1974.10 
NRGLEV= 4 
NTABLS* 3 
CUTOFF- 23.99 
I SKIP- 1 
KPRINT= 1 


T-2 

Parameter Output 






protons (E> 5.00 MEV) **4t** 


LONG* 

L AT. 

ALT. 

FIELD 

LINE 

D^BIT 

POSITIONAL 

TIME-INTEG 




-B- 

-L- 

TIME 

FLUX 


PSTNL FLUX 

(DEG) 

(DEG) 

(KM) 

(GAUSS) 

( E • R • > 

( HRS) 

#/CM**2/SEC 



- 99 • 7 3 

0 *02 

1994 *0 

0*141 10 

1.36 

0. 0 





- 1 00 • 2 3 

5*70 

1994*2 

0*15100 

1.42 

0. 03333 

2. 339E 

04 

2.806E 

06 

-1 00* 75 

1 1 *37 

1994.7 

0.16342 

1.51 

0. 06667 

1 . 69 2E 

04 

2 . 03 IE 

06 

-101*23 

17*05 

1995*5 

0.17749 

1.65 

0. 10000 

1. 01 IE 

04 

1 .213E 

06 

- 1 01 *73 

22*72 

1996*6 

0.19234 

1.64 

0.13333 

4.953E 

03 

5.944E 

05 

- 1 02 • 2 4 

28.3 9 

1997.9 

0.20714 

2.12 

0.16667 

3.589E 

03 

4.307E 

05 

-102*74 

34*06 

1999.5 

0.22114 

2.51 

0.20000 

9.497E 

02 

1 • i 40E 

05 

-103*24 

39*73 

2001 • 1 

0.23371 

3.09 

0. 23333 

5.265E 

01 

6. 318E 

03 

-103. 74 

45.40 

2002 *8 

0.24437 

3.96 

0.26667 

0*0 


0. 0 


-104*24 

51 *06 

2004 .4 

0*25282 

5.36 

0.30000 

0.0 


0*0 


-104.74 

56.73 

2006*0 

0*25896 

7. 78 

0.35333 

0. 0 


0.0 


-105*24 

62.39 

2007 .3 

0*26294 

1 2. 43 

0.36667 

0.0 


0* 0 


-105* 74 

68*05 

2008 *5 

-1.00000 

- 1 . 00 

0.40000 

0.0 


0 *0 


- 1 06 . 2 5 

73 . 7 1 

2009*4 

-1 .00000 

-1 .00 

0.43333 

0. 0 


0.0 


-106* 75 

79 .37 

2010*0 

-1.00000 

- 1 • 00 

0.46667 

0.0 


0.0 


-107*24 

85.03 

2010 .2 

-1.00000 

-1 .00 

0.50000 

0.0 


0.0 


72. 1 9 

89*31 

2010.1 

-1*00000 

-1.00 

0.53333 

0.0 


0*0 


71 .74 

83*65 

2009 .6 

-1.00000 

-1 . 00 

0.56567 

0.0 


0.0 


71 *24 

77.99 

2008.8 

0.26230 

1 2. 01 

0. 60000 

0. 0 


0.0 


70*74 

72 .33 

2007.6 

0*25992 

7.58 

0.63333 

0. 0 


0. 0 


70.24 

66.67 

2006*2 

0.25619 

5.25 

0.66667 

0.0 


0.0 


69*74 

61 *00 

2004*5 

0.250 76 

3.87 

0. 70000 

0.0 


0*0 


69* 24 

55*34 

2002.6 

0.24343 

3.01 

0. 73333 

8. 072E 

01 

9.686E 

03 

68.74 

49.67 

2000 .6 

0.23417 

2.43 

0. 76667 

7.752E 

02 

9.303E 

04 

68*24 

44.00 

1998.6 

0.22314 

2. C4 

0.80000 

1.557E 

03 

1 . 868E 

05 

67.74 

38 *32 

1996.6 

0.21069 

1 .76 

0.83333 

1.659E 

03 

1 .990E 

05 

67. 23 

32.65 

1994.7 

0.19741 

1.57 

0.86667 

4. 108E 

03 

4.930E 

05 

66.73 

26*97 

1993 .0 

0*18415 

1.43 

0.90000 

6. 102E 

03 

7.322E 

05 

66.23 

21 .29 

1991 .5 

0.17203 

1 .34 

0.93333 

8.266E 

03 

9.919E 

05 

65.7 3 

15.61 

1990 .3 

0.16230 

1.29 

0. 96667 

1.031E 

04 

1 . 237E 

06 

65*23 

9.93 

1989.4 

0.15608 

1 .26 

1 .00000 

1 . 1 5 5E 

04 

1 . 386E 

06 

64.7 3 

4.25 

1986.9 

0.15392 

1 .27 

1 .03333 

1.35 0E 

04 

1 • 62 OE 

06 

64* 23 

-1 *44 

1986.8 

0.15554 

1 .31 

1 . 06 6 6 7 

1 . 73 OE 

04 

2. 076E 

06 

63.72 

-7*12 

1989.1 

0.16002 

1.38 

1 .10000 

1.51 3E 

04 

1 .816E 

06 

63.22 

-12*80 

1989.8 

0.16619 

1.49 

1 • 13333 

1.491E 

04 

1 .789E 

06 

62*72 

-13*48 

1990.9 

0.17308 

1 .65 

1 • 16667 

1.205E 

04 

1.4 46E 

06 

62. 22 

-24*16 

1992.2 

0.18010 

1.87 

1.20000 

8. 208E 

03 

9.850E 

05 

61 .72 

-29.84 

1993 .8 

0.18701 

2.17 

1.23333 

7.568E 

03 

9. 081E 

05 

61 *22 

-35*52 

1995-7 

0.19380 

2.57 

1 .26667 

1.777E 

03 

2. 133E 

05 

60.72 

-41.19 

1997.6 

0.20057 

3. 14 

1.30000 

7.793E 

01 

9.352E 

03 

60*22 

-46.86 

1999.7 

0.20739 

3.94 

1.33333 

0. 0 


0. 0 


59.71 

-52*53 

2001 .7 

0.21427 

5.10 

1.36666 

0.0 


0.0 


59 .21 

-58.20 

2003*6 

0.22116 

6.67 

1.40000 

0.0 


0.0 


58*71 

-63*86 

2005.4 

0.22789 

9*61 

1.43333 

0.0 


0. 0 


58*21 

-69*53 

2007.0 

0.23425 

1 3. 83 

1.46667 

0.0 


0.0 


57.71 

-75.19 

2008.3 

-1 .00000 

-1 .00 

1.50000 

0.0 


0.0 


57.21 

-80.85 

2009.3 

-1 .00000 

-1.00 

1.53335 

0.0 


0.0 


56.71 

-86.51 

2009*9 

-1.00000 

-1,00 

1.56667 

0.0 


0.0 


-125*80 

-87*83 

2010 *2 

-1.00000 

-1 .00 

1.60000 

0.0 


0.0 


-124*30 

-82.17 

2010*2 

-1.00000 

-1 .00 

1.63333 

0.0 


0.0 


-124.80 

-76.51 

2009.7 

0.247E5 

8.25 

1.66667 

0. 0 


0. 0 


- 1 25 . 3 0 

-70*85 

2009.0 

0.24372 

5*93 

1. 70000 

0. 0 


0. 0 


-125*60 

-65*19 

2007 .9 

0.23761 

4.45 

1.73333 

0. 0 


0. 0 


- 1 26 • 30 

-59*53 

2006*6 

0.22970 

3.48 

1 . 76667 

5.619E 

00 

6. 743E 

02 

-126*80 

-53.87 

2005*2 

0.22024 

2.85 

1 .80000 

3.318E 

02 

3.981E 

04 

-127*31 

-48.20 

2003.5 

0.20956 

2.37 

1 .83333 

2.961E 

03 

3. 553E 

05 

-127.81 

-42 *53 

2001 *9 

0*19302 

2.04 

1 *86666 

5.884E 

03 

7-.-06^0E 

05 

- 1 26 • 3 1 

-36*87 

2000*2 

0.18603 

1 .80 

1.90000 

6.74 9E 

03 

8. 099E 

05 

-128*81 

-31 .20 

1998*6 

0.17411 

1 .62 

1.93333 

1.168E 

04 

1 .401E 

06 

-129.31 

-25.52 

1997*2 

0.16288 

1.50 

1 .96667 

1.700E 

04 

2. 040E 

06 

-129.81 

-1 9*85 

1996*0 

0.15309 

1 .41 

2.00000 

2.143E 

04 

2.572E 

06 

-130.31 

-14.18 

1995 .0 

0.14558 

1 .35 

2. 03333 

2.601E 

04 

3 . 1 21 E 

06 


ORBITAL 

FLUX(SUM) 


2.806E 

06 

4.837E 

06 

6.0SOE 

06 

6.644E 

06 

7.075E 

06 

7.189E 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.19SE 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.195E 

06 

7.205E 

06 

7.298E 

06 

7.485E 

06 

7.684E 

06 

8.177E 

06 

8.909E 

06 

9.901E 

06 

1 . 1 1 4E 

07 

1.2S2E 

07 

1.414E 

07 

1 .622E 

07 

1.804E 

07 

1.983E 

07 

2.127E 

07 

2.226E 

07 

2.316E 

07 

2.338E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.339E 

07 

2.343E 

07 

2.378E 

07 

2.449E 

07 

2.530E 

07 

2.670E 

07 

2.874E 

07 

3.131E 

07 

3.443E 

07 


T-3 

Beginning of Running Printout 



-84. 17 

-65.24 

2008.1 

0.20507 

3.42 

22.89999 

1.952E 

01 

2.343E 

03 

3.400E 

08 

-84.67 

-59.5 8 

2006*9 

0.19271 

2.76 

22.93330 

7.61 3E 

02 

9. 136E 

04 

8.401E 

08 

-85 • 1 7 

-53 .91 

2005.4 

0.18030 

2.31 

22. 96666 

6. 676E 

03 

8.012E 

05 

8.409E 

08 

-85.68 

-48.25 

2003.8 

0.16828 

2.00 

22.99998 

1 . 58 OE 

04 

X «896E 

06 

8.42aE 

08 

-86.18 

-42.56 

2002.2 

0.15704 

1.78 

23.03331 

2.085E 

04 

2.502E 

06 

8.453E 

03 

-86*68 

-36.91 

2000 .5 

0.14696 

1.62 

23.06667 

3.203E 

04 

3.844E 

06 

8.491E 

08 

- 87 . 1 8 

-31 .25 

1999.0 

0.13844 

1 .51 

23.09999 

4.250E 

04 

5.100E 

06 

8.542E 

08 

-87.68 

-25.5 7 

1997.5 

0.13194 

1.43 

23. 13332 

4.992E 

04 

5.991E 

06 

8.602E 

08 

-88. 18 

-19 .90 

1996.3 

0.12800 

1.38 

23.16664 

5. 73 2E 

04 

6.878E 

06 

8.671E 

08 

-88.68 

-14.23 

1995.4 

0.12716 

1.35 

23.20000 

5.855E 

04 

7. 026E 

06 

8.741E 

08 

-89.18 

-8 .55 

1994.7 

0.12983 

1 .35 

23.23331 

5. 195E 

04 

6.234E 

06 

8.803E 

08 

-89*69 

-2.88 

1994.3 

0.13610 

1.37 

23.26666 

3.948E 

04 

4. 738E 

06 

8.851E 

08 

-90.19 

2.80 

1994.3 

0.14567 

1.42 

23.29999 

2.863E 

04 

3. 436E 

06 

' 8.885E 

08 

-90.69 

8.47 

1994.7 

0.15790 

1.50 

23. 33331 

2.068E 

04 

2. 482E 

06 

8.910E 

08 

-91.19 

14.15 

1995.3 

0.17196 

1 .62 

23.36664 

1.272E 

04 

l.526E 

06 

8.925E 

08 

-91.69 

19.82 

1996.3 

0.18700 

1 .80 

23.39999 

6. 474E 

03 

7 . 769E 

05 

8.933E 

06 

-92.19 

25.49 

1997.5 

0.20217 

2.05 

23.43330 

4.860E 

03 

5.833E 

05 

8.939E 

08 

-92.69 

31 .16 

1998.9 

0.21671 

2*41 

23.46666 

2.43SE 

03 

2.922E 

05 

8.942E 

08 

-93.19 

36.83 

2000.4 

0.22995 

2.92 

23.49998 

1.340E 

02 

1 .608E 

04 

8.942E 

08 

-93.70 

42 .50 

2002.1 

0.24136 

3.71 

23.53331 

0. 0 


0. 0 


8.942E 

08 

-94.20 

48.1 7 

2003.7 

0.25058 

4.94 

23.56667 

0.0 


0.0 


8.942E 

08 

-94.70 

53.83 

2005.3 

0.25745 

7.02 

23.59999 

0.0 


0.0 


8.942E 

08 

-95.20 

59.50 

2006.8 

0.26205 

1 0.93 

23.63332 

0.0 


0.0 


8.942E 

08 

-95.70 

65.16 

2008.0 

0.26466 

19.48 

23.66664 

0. 0 


0.0 


8.942E 

08 

-96.20 

70.82 

2009.0 

-1.00000 

-1 .00 

23. 70000 

0.0 


0.0 


8.942E 

08 

-96. 70 

76 .46 

2009.6 

-1.00000 

- 1 . 00 

23. 73331 

0.0 


0.0 


a.942£ 

08 

-97.20 

82.14 

2010.1 

-1.00000 

-1. 00 

23. 76666 

0.0 


0.0 


8.942E 

08 

-97.70 

8 7.80 

2010.2 

-1.00000 

-1 .00 

23.79999 

0.0 


0.0 


8.942E 

08 

81 .78 

86.54 

2009.9 

-1.00000 

-1 .00 

23.83331 

0.0. 


0.0 


8.942E 

08 

81 .29 

80 .88 

2009.2 

0.26379 

1 5.41 

23. 86664 

0.0 


0.0 


8.942E 

08 

80.79 

75.22 

2006.1 

0.26251 

9.18 

23.89999 

0.0 


0.0 


8.942E 

08 

60.29 

69 .56 

2006.8 

0.26016 

6.12 

23.93330 

0.0 


0.0 


8.942E 

06 

79.78 

63 .90 

2005.2 

0.25625 

4.40 

23.96666 

0.0 


0.0 


8.942E 

08 


percent of total lifetime spent inside and outside trapped particle radiation belt 

INNER ZONE (1*0 <= L < 2*6) : 51*60 % 

OUTSIDE TRAPPING REGION ( I . 0 <= L < 1 . 1 > : 0.0 % 

INSIDE TRAPPING REGION (1.1 <= L < 2.8) : 51.60 % 

OUTER ZONE <2.8 <= L <= 11.0) : 26.84 % 

EXTERNAL (L > II.O) : 2 1.56 % 

TOTAL ORBIT TIME IS : 24.00 HOURS 


T-4 

Final Part of Running Printout with Percent Time Table 



* Jiafc 4c ***:i^*Jft ****** ******* *:k ** Ic ti *********** ******** **** *************************************** ********************************** 

* SOFIP : SHORT ORBITAL FLUX INTEGR. PRCGRAM FOR STANDARD NSSDC PROTON AND ELECTRON ENVIR. MODELS ISPECIES CONSIDERED SEPARATELY! ♦ 

* MAGNETIC PARAMETERS B AND L COMPUTED tolTH GEOMAGN. FIELD MODEL 5: IGRF 1965.0 80-TERM 10/68 ♦ COEFF. UPDATED TO: 1974.1 * 

* PROJECT :TEST 90/2000 * INCLIN= 90CEG ♦ PERIG= 2000KM ♦ AP0G= 2000KM ♦ B/L TAPE=TD5116 * PERIOD= 2.120HRS * SOLA^ MINIMUM * 

* FDR INFORMATION OR EXPLANATION CONTACT E. G. ST ASSI NO^OUL OS AT NASA- GSFC .CODE 601, GREEN8ELT, MARYLAND 20771, TEL .( 30 1 ) -344-8067 * 

* **it*k-*******r^r9i******'jit*K:*^**:^* ******** * 4i******** ************************************ ** ********************************************* 

**************** PROTCNS **************** 

******* ************A ******************* ********** 


*************** CCMPCSITE orbit SPECTRUM *************** 

ENERGY AVERAGED DIFFERENCE AVERAGED DIFFE- 

LEVELS INTEGRAL FLUX INTEGRAL FLUX RENTTAL FLUX 

>(MEV> iO</CM+*2/SEC #/CM**2/SEC/DE #/CM»*2/SEC/KFV 

2.000 1.646E 04 2.454E 03 0.0 

2.000 1.401E 04 1.667E 03 0.0 

4.000 1.234E 04 1.218E 03 0.0 

5.000 1.112E 04 9.428E 02 0.0 

6.000 1.018E 04 1.615E 03 0.0 

e.OOO 8.564E 05 1.247E 03 0.0 

10.00 7.318E 03 1.642E 03 0.0 

15.00 5.675E 03 1.203E 03 0.0 

20.00 4.472E 03 4.420E 02 0.0 

25.0 0 4.0 30E 03 3.913E 0 2 0.0 

30.00 2.639E 03 2.295E 02 0.0 

25.00 3.409E 03 2.138E 02 0.0 

40.00 3.195E 03 1.993E 02 O.Q 

45.00 2.996E 03 1.360E 02 0.0 

50.00 2.810E 03 1.528E 02 0.0 

55.00 2.657E 03 1.441E 02 0.0 

60.00 2.513E 03 2.642E 02 0.0 

70.00 2.249E 03 2.35SE 02 0.0 

60.00 2.014E 03 2.100E 02 0.0 

90.00 1.804E 03 1.874E 02 0.0 

100.0 1.616E 03 3.818E 02 0.0 

125.0 1.234E 03 2.901E 02 0.0 

150.0 9.444E 02 2.209E 02 0.0 

175.0 7.235E 02 1 . 68SE 02 0.0 

200.0 5.549E 02 2.143E 02 0.0 

250.0 3.407E 02 1.311E 02 0.0 

200.0 2.096E 02 8.042E 01 0.0 

550.0 1.292E 02 4.944E 01 0.0 

400.0 7.973E 01 4.925E 01 0.0 

500.0 3.047E 01 3.047E 01 0.0 
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Output Table 1 without Differential Spectrum, Solar Protons, or Exposure Index - Protons 



4 e ♦♦ *♦ >ft5jc**<c ** *** * 4 : *«*5je*#*)iK#*’*'«:***********^**=<‘************’^'**=<^*** ’*'*** ****** 

* SQFIP : SHORT ORBITAL FLUX INTEGR. PROGRAM FOR STANDARD NSSOC PROTON AND ELEC TRUN ENV IR« MODELS ( SPEC IE S CONS I DERED SEP AK AT EL Y ) * 

* MAGNETIC PARAMETERS 0 AND L COMPUTED KlITH GEOMAGN* FIELD MODEL 5: IGRF 1965*0 80—TERM 10/63 * COEFF* UPDATED TO. 1974*1 * 

t PRo!jicT sloFlP TEST * INCLIN= 900ES ♦ PER1G= 2000KM * APOG= 2000KM » B/L TAPE=TD5116 * PER ICD= 2ol20HRS »■ SOL AR M N IMUM * 

* INFORMATION OR EXPLANATION CONTACT E*G* STASS INOPOULOS AT N ASA-GSFC » CODE 601* GREENBELT* MARYLAND ii0771* TEL •< 30 1 I -344-806 7 * 

:ti^*^i****^*^* ***^^4 ^*t¥*^******** **♦*♦♦♦*♦♦»***=»********************************************************************** 


ELECTRONS LO ❖sJ' *********** * 

******************************************* ****** 


*************** CCMFOSITE OKdl T SPECTRUM *************** 

ENERGY AVERAGED DIFFERENCE AVERAGED DIFFE- 

LEVELS INTEGRAL FLUX INTEGRAL FLUX RENT lAL FLUX 

X MEV) #/CM**2/SEC #/CM**2/SEC/ DL #/CM * *2/SEC/ K£ V 


• 1000 

7.454E 

06 

4. 716E 

06 

0.0 

• 2000 

2.73SE 

06 

1.521E 

06 

0.0 

• 3000 

1. 217E 

06 

5.C30E 

05 

0.0 

• 4000 

6.538E 

05 

2. 962E 

05 

0.0 

.5000 

3.576E 

05 

1 . 1 1 1 E 

05 

0.0 

• 6000 

2. 465E 

C5 

7.425E 

04 

0.0 

• 7000 

1. 722E 

05 

4p388E 

04 

0 o 0 

.8000 

1. 284E 

05 

2.718E 

04 

0.0 

.9000 

1.012E 

C5 

2. 1 17E 

04 

0.0 

1*000 

8.001E 

04 

3.057E 

04 

0«0 

1.250 

4.944E 

04 

1 .677E 

04 

0.0 

l.SOO 

3. 067E 

04 

1.046E 

04 

0.0 

1.750 

2. 021E 

C4 

6.862E 

03 

Oo 0 

2.000 

1.334E 

04 

3. 930E 

03 

0.0 

2.250 

9.414E 

03 

2.740E 

03 

0.0 

2.500 

6.674E 

C3 

3o037E 

03 

Oo 0 

2.750 

3.637E 

C3 

1.449E 

03 

0.0 

3.000 

2. 188E 

03 

1.003E 

03 

0.0 

3.250 

1 . 185E 

C3 

4. 794E 

02 

OoO 

3.500 

7.055E 

02 

2.620E 

02 

0.0 

3.750 

4.435E 

02 

1. 567E 

02 

0.0 

4 .000 

2.868E 

02 

1.484E 

02 

0.0 

4.250 

1. 384E 

02 

7. 141E 

01 

0.0 

4.500 

6.702E 

01 

3.451E 

01 

0.0 

4.750 

3.252E 

01 

1.669E 

01 

0.0 

5.000 

I.583E 

01 

io426E 

01 

0.0 

5.500 

1 .550E 

00 

1.525E 

00 

0.0 

6.000 

2. 509E 

-02 

2.509E-02 

0.0 

6.500 

0. 0 


Oo 0 


0.0 

7.000 

0. 0 


0.0 


0. 0 
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, or Exposure Index - 


Output Table 1 without Differential Spectrum, Solar Protons 


Electrons 



46 SOFI 

♦ PRO JECT*" : TESt'^90/20 00 *'^*INClTn= 90DEG * PERIG= 200okM ♦ APOG= 2000KM # B/L. TAPE=fosil6 # PERIOC= 2.120HRS * SOLAR MAXIMUM 

FOR INFORMATION DR EXPLANATION CONTACT E«G« ST ASS I NO ^OUL 3S AT NASA-GSFC tCOOE 601# 6REENBELT# MARYL^ND 20771 # ^ . f 

it 41*45 ♦♦♦♦♦♦ + *41 * ** **+ ♦♦♦♦♦*♦♦***♦***♦** 


A* 41 54 c* 4 cA 41 **41***41*41 *41*** 41 


**A****A******* COMPOSITE ORBIT SPECTRUM *************** 


ENERGY AVERAGED DIFFERENCE 

LEVELS INTEGRAL FLUX INTEGRAL FLUX 
XMEV) #/CM**2/SEC #/C M ** 2/S EC/D E 


AVERAGED DIFF- 
RENTI AL FLUX 
#/CM**2/SEC/KEV 


2*000 

1 .5042 

04 

2. 134E 

0 3 

2.242E 00 

3.00 0 

1 .291E 

04 

1 .468E 

03 

1.7 2 3E 00 

4.000 

i .144E 

04 

1 .089E 

03 

1.297E 00 

5. 000 

1 .035E 

04 

8. 553E 

02 

1.0 0 2E 00 

6. 000 

9.494E 

03 

1 .536E 

03 

8.395E-01 

8.000 

7.958E 

03 

1 . 190E 

03 

6.46 lE-Oi 

10.00 

6.768E 

03 

1 .519E 

03 

4.689E-01 

15.0 0 

5.249E 

03 

1 . 1 1 3E 

03 

2.529F-01 

20.0 0 

4. 136E 

03 

4. 300E 

02 

1 .471E-01 

25. 0 0 

3.706E 

03 

3. 791E 

02 

9.040E-02 

30.0 0 

3.327E 

03 

2. 165E 

02 

5.933E-02 

35.00 

3.1 1 oe 

03 

2.013E 

02 

4.S29E-02 

40.0 0 

2.909E 

03 

1 . 673E 

02 

3.901E-02 

45 .0 0 

2.722E 

03 

1.7 44E 

02 

3.52 2E-02 

50.0 0 

2.547E 

03 

1 .386E 

02 

3.11 7E-02 

55.00 

2.409E 

03 

1. 307E 

02 

2.762E-02 

60.0 0 

2.278E 

03 

2.39 7E 

02 

2. 53 5E-02 

70.00 

2.038E 

03 

2. 136E 

02 

2.259E-02 

80.00 

1 .825E 

03 

1 .905E 

02 

2.016E-02 

90.0 0 

1.634E 

03 

1 • 700E 

02 

1.79 6E-02 

100.0 

1 .464E 

03 

3. 468E 

02 

1.598E-02 

125.0 

1 .1 17E 

03 

2. €33E 

02 

1.20 5E-02 

150.0 

8.540E 

02 

2.003E 

02 

9. 15 6E-03 

175.0 

6.537E 

02 

1.527E 

02 

6.906E-03 

200.0 

5.010E 

02 

1 .934E 

02 

5.146E-03 

250. 0 

3.076E 

02 

1 • 1 83E 

02 

3.02 0E-03 

300.0 

1 .893E 

02 

7.255E 

01 

1. 83 3E- 03 

350.0 

1 .1672 

02 

4. 461E 

01 

1. 12 7E-03 

400.0 

7.213E 

01 

4.448E 

Cl 

6.939E-04 

500 • 0 

2.765E 

0 1 

2.765E 

01 

2.64 9E-G4 


**** SOLAR PROTONS **** 


FOR TAU=12#O=90: NALE=l 
WITH GEQMAG SHIELDING 
lEXPOSR FACTQR=0*34) 

ENERGY TOTAL 

LEVELS FLUENCE 

XMEV) #/CM**2 


*♦ EXPOSURE index: ENERGY>5«000 MEV ** 


I NTENSITY 
RANGES 
#/CM**2/SEC 


EXPOSURE TOTAL » OF 
DURATION ACCUMULATED 
(HOURS) PARTICLES 


1 0. 

S.741E 

09 

ZERO FLUX 

10.100 

0. 0 


20. 

3.937E 

09 

1 

.EO-1 .El 

0.100 

2.749E 

03 

30. 

2.699E 

09 

1 

.El-1 .E2 

1 .000 

1.686E 

05 

40. 

1 .85 IE 

09 

1 

.E2-1.E3 

1 .333 

2.22SE 

06 

50. 

1 .269E 

09 

1 

.E3-1 .E4 

5.000 

a. 745E 

07 

60. 

8. 70 IE 

08 

1 

.E4-1 .E5 

6.267 

7.407E 

06 

70. 

5.966E 

08 

1 

.E5-1 .E6 

0.167 

6.368E 

07 

80. 

4.091E 

08 

1 

.E6-1 .E7 

0*0 

0.0 


90* 

2.805E 

08 

1 

.E7-0VER 

0.0 

0.0 


100. 

1.923E 

08 


TOTAL 

23.964 

8.942E 

08 

110. 

1.319E 

08 






120. 

9.042E 

07 






1 30. 

6.200E 

07 






140. 

4.251E 

07 






1 50. 

2.915E 

07 






1 60. 

1 .999E 

07 






1 70. 

1 .370E 

07 






180. 

9.396E 

06 






190. 

6.443E 

06 






200. 

4.41 8E 

06 
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Output Table 1 with Differential Spectrum, Solar Protons, and Exposure Index ~ Protons 
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♦ SOFIP : SHORT ORBITAL FLUX INTEGR, PROGRAM FOR STANDARD NSSDC PROTON AND iiLECTRON ENVIR, MODELS (SPECIES CONSIDERED SEPARATELY) ♦ 

♦ MAGNETIC PARAMETERS B AND L COMPUTED WITH GEOMAGN* FIELD MODEL S: IGRF 1965.0 80-TERW 10/68 ♦ COEFF* UPDATED TQ: 1974.1 ♦ 

* PROJECT :SOFIP TEST ♦ INCLIN= 900EG ♦ PERIG= 2000KM ♦ AP0G= 2000KM * B/L TAPE=TD5ll6 * PER ICD= 2®120HRS * SOLAR MAXIMUM ♦ 

* FDR INFORMATION OR EXPLANATION CONTACT E.G. STASS INOPOULOS AT NASA- GSFC * CODE 601, GREENBELT, MARYLAND 20771, TEL . ( 30 1 )-344-8067 ^ 

*♦ ♦♦♦♦*♦ *4 ♦♦♦★♦:(J***** 

ELECTRONS Hi #*+******jjc****4t* 
Hiitciiitifi^^iiii*nf***m********^^****^**********^********** 


*♦♦♦ SOLAR PROTONS *♦♦♦ 

FOR TAU-12,Q=90: NALE=1 

COMPOSITE ORBIT SPECTRUM ♦:»*♦♦♦♦♦♦♦♦♦*♦* WITH GE0MA6 SHIELDING ** EXPOSURE INDEX: ENERGY>®5000 MEV *'* 

(EXPOSR FACTl)R=0.34) 


ENERGY 

AVERAGED 

DIFFERENCE 

AVERAGED 

1 DIFF- 

ENERGY 

TOTAL 


INTENSITY 

EXPOSURE 

TOTAL » 

OF 

LEVELS 

INTEGRAL 

FLUX 

INTEGRAL 

FLLX 

RENTIAL 

FLUX 

LEVELS 

FLUENCE 


RANGES 

DURATION 

ACCUMULATED 

XMEV) 

#/ CM ♦♦2/ SEC 

«/CM** 2/SEC/DE 

tf/CM**2/SEC/KEV 

XMEV) 

#/CM**2 


tf/CM*+2/SEC 

(HOURS) 

PARTI CLES 

• 1000 

1.511E 

07 

8.096E 

06 

1. 175E 

05 

10. 

5.74 IE 

09 


ZERO FLUX 

E.333 

0. 0 


• 2000 

7.012E 

06 

4.037E 

06 

S.685E 

04 

20o 

3.937E 

09 


1 .EO-l.El 

0.033 

2. 1 71E 

02 

• 3000 

2.9 76E 

06 

1.821E 

06 

2.557E 

04 

30. 

2.699E 

09 


1. El-1, E2 

0.200 

1.575F 

04 

• 4000 

1 . 1 54E 

06 

6.924E 

05 

9.787E 

03 

40. 

1.051E 

09 


1 .E2-1.E3 

0.333 

5.381E 

05 

.5000 

4.62IE 

05 

1.672E 

05 

3.515E 

03 

50. 

1.269E 

09 


1 .E3-1.E4 

1.733 

2.886E 

07 

t.6000 

2.949E 

05 

1.014E 

05 

1.461E 

03 

60. 

6.70i£ 

OS 


1 .E4-1.E5 

3.267 

5.272E 

08 

• 7000 

1.934E 

05 

5.586E 

04 

7.461E 

02 

70. 

5.966E 

08 


I.E5-1.E6 

9.367 

1.317E 

10 

• 8000 

1.376E 

05 

3.302E 

04 

4.435E 

02 

80. 

4.091E 

03 


lo E6-1.E7 

3.700 

2.619E 

10 

• 9000 

1 .046E 

C5 

2.455E 

04 

2.844E 

02 

90 • 

2.805E 

08 


1. £7- OVER 

0.0 

0.0 


1*000 

8.001E 

04 

2.8 59E 

04 

1 . 869E 

02 

100. 

1.923E 

08 


TOT AL 

23.96 4 

3.993E 

10 

1.250 

5.143E 

04 

1.799E 

04 

9.252E 

01 

110. 

1.31 9E 

08 






1.500 

3.343E 

04 

1.0 33E 

04 

5.388E 

01 

120. 

9.042E 

07 






1.750 

2.311E 

04 

7.046E 

03 

3.426E 

01 

130. 

6.200E 

07 






2*000 

1.606E 

04 

4.479E 

03 

2.256E 

01 

140, 

4.251E 

07 






2.250 

1.158E 

C4 

3.1 85E 

03 

1.673E 

01 

150. 

2.915E 

07 






2.500 

8.395E 

03 

3.395E 

03 

1.303E 

01 

160. 

1.999E 

07 






2.750 

5.000E 

03 

1.735E 

03 

8.791E 

00 

170. 

1.370E 

07 






3.000 

3.265E 

03 

9.799E 

02 

5. 154E 

00 

180. 

9.396E 

06 






3«250 

2.285E 

03 

5.342E 

02 

2.950E 

00 

190. 

6.443E 

06 






3.500 

1.751E 

03 

3.551E 

02 

1.833E 

00 

200. 

4.418E 

06 






3.750 

1.396E 

03 

2.645E 

02 

1.395E 

00 









4.000 

1. 131E 

03 

3.415E 

02 

1.256E 

00 









4.250 

7.897E 

02 

2.3 70E 

02 

1.059E 

00 









4.500 

5.527E 

02 

1 .6 53E 

02 

7.869E- 

01 









4.750 

3.874E 

02 

1 • 1 55E 

02 

6. 0046- 

01 









5. 000 

2.7I9E 

02 

1.9 17E 

02 

4.885E- 

01 









5.500 

8.018E 

01 

5.648E 

01 

2. 194E- 

01 









6. 000 

2.370E 

01 

2.159E 

01 

6.402E- 

02 









6.500 

2.1 13E 

00 

2.1 13E 

00 

9.E18E- 

03 









7.000 

0.0 


0.0 


0.0 
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Output Table 1 with Differential Spectrum, Solar Protons, and Exposure Index - Electrons 



* SOFIP : SHORT ORBITAL FLUX INTE3R-. PROGRAM FOR STANDARD NSSDC PROTON AND ELECTRON ENVIR« MODELS ( SPECIES CONSIDERED SEPARATELY} * 

* MAGNETIC PARAMETERS B AND L COMPUTED WITH GEOMAGN, FIELD MODEL SI IGRF 1965.0 80-TERM 10/68 ♦ COEFF. UPDATED TO: 1974*1 * 

* PROJECT ITEST 90/2000 ♦ INCLIN= 90DEG * PERIG= 2000KM if APOG= 2000KM ♦ B/L TAPE=T05116 » PERIOD= 2.120HRS * SOLAR MAXIMUM ♦ 

* FOR INFORMATION OR EXPLANATION CONTACT E.G* STASSI NO^OUL 3S AT NASA-GSFC .CODE 601. GREENBELT. MARYLAND 20771. TE^ .( 30 11 -344-8067 4 

4 ( 4(4 49 ^ 44 4444(4 44 444 44 44444444444444444 4 4444 4 44444 44 4 44 444 4 4 444 4 4 4 4444444 4 4444 4 44 444 444 44 44 44444444 4 

444444444444444444444444 PROTCNS 444444444444444444444444444 

44 TABLE OF PEAK AND TCTAL FLUXES PER PERIOD : ENERGY > 5*0 MEV 44 

444444444444444*444444444444444444444444444 4444444444444444444444444 


PERIOD 

PEAK FLUX 

PCSITICN 

AT WHICH ENCOUNTERED 

ORBIT TIME 

FIELO(B) 

LINE(L) 

TOTAL 1 

FLUX 

NUMBER 

ENCOUNTERED 
4 /CM* 4 2/SEC 

LCNGI TUDE 
(DEG) 

LATI TUDE 
( DEG) 

ALTITUDE 
( KM) 

(HOURS) 

(GAUSS) 

(E.R. ) 

PER ORBIT 
4/CM 4* 2/OR BIT 

1 

a.253E 

04 

- 131*316 

-2.83 

1 993.99 

2.10000 

0. 140 41 

1 .31 

4.552E 

07 

2 

3*599E 

04 

31*637 

-9.95 

1989.39 

3.23333 

0.14089 

1 .46 

6. 160E 

07 

3 

6.541E 

04 

-0.451 

-12.79 

1989.68 

5.36666 

0.12769 

1 .50 

8.547E 

07 

4 

1 .059E 

05 

-33.040 

-21.31 

1991 .28 

7.53333 

0.11669 

1 *50 

1.200E 

08 

S 

9*788E 

04 

-64.626 

-18.46 

1990.65 

9.63333 

0.11720 

1 *39 

9.174E 

07 

6 

4.973E 

04 

-96*212 

-15.62 

1990. 13 

11.73333 

0* 13084 

1 .35 

5.599E 

07 

7 

3 •321E 

04 

-127.298 

-7.08 

1 988.98 

13.80000 

0.13994 

1 *31 

6.099E 

07 

8 

5.778E 

04 

6.577 

-11.38 

1994. 78 

16*86664 

0.13019 

1 .49 

7.891E 

07 

9 

9.856E 

04 

-25.511 

-8.54 

1994.54 

16.99998 

0.11706 

1 .40 

1 .ia7£ 

08 

10 

1 *099E 

05 

-56.596 

-17. 06 

1 995. 76 

21.06667 

0*11493 

1 *40 

1 « 069E 

08 

11 

S.aSSE 

04 

-88.683 

-14.23 

1995.36 

23.20000 

0.12716 

1.35 

6.306E 

07 
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Output Table 2 


Protons 



♦ SDFIP S SHORT ORBITAL FLUX INTEGR. PROGRAM FOR STANDARD NSSDG PROTON AND ELECTRON ENVIR. MODELS ( SPEC iE S CONSl DERED SEPARATELY) 

♦ MAGNETIC PARAMETERS 8 AND L COMPUTED WITH GEOMAGN* FIELD MODEL 5: IGHF 1965*0 SO-TERM 10/68 * COEFF* UPDATED TO: 1974.1 

♦ PROJECT :SOFIP TEST ♦ INCLIN= 900EG ♦ PERI G^= 2 000KM * APOG= 2000KM *■ b/L TAPC=TD5116 * PER IOD= 2o 120HRS * SOLAR MAXIMUM 

♦ FOR INFORMATION OR EXPLANATION CONTACT E*G* STASS TNOPGULOS AT NASA-GSFC , CODE 601* GREENBLLT* MARYLAND 20771, T EL . ( 30 1 ) -344 -8067 

:if:iHtiitt4********Mf**^ ******* ELECTRONS HI ********^:t^***************** 

** TABLE OF PEAK AND TOTAL FLUXES PER PERIOD : ENERGY > 0.50 MEV ** 

itt4c*nc**^** *************************** ***im*^ ************ ********** 


PER lOD 

PEAK FLUX 

POSITION 

AT WHICH ENCOUNTERED 

ORBIT TIME 

F IELD{ 8) 

LINE(L ) 

TOTAL 

FLUX 

NUMBER 

ENCOUNTERED 
il^/CM* ♦2/SEC 

LONGITUDE 

(DEG) 

LATITUDE 

(DEG) 

ALT (TUDE 
(KM) 

(HOURS) 

( GAUSS ) 

(E.R. ) 

PER ORBIT 
«/CM4*2/0ReiT 

1 

1 .980E 

06 

-129.812 

-19.85 

1995.97 

2.00000 

0.15309 

1.41 

2.e«9E 

09 

2 

2.216E 

06 

31*637 

-9*95 

1939.39 

3.23333 

0. 14089 

1. 46 

3.260E 

09 

3 

2.940E 

06 

0*051 

-7.11 

1 989.01 

5.33333 

0.12639 

1.42 

3.779E 

09 

4 

3.536E 

06 

-32.538 

- 15.63 

1 990. 08 

7.50000 

0.11529 

1.44 

4.472E 

09 

5 

3.309E 

06 

-65.127 

-24.14 

1 992.0 1 

9. 6666 7 

0.11888 

1.43 

3o650E 

09 

6 

2.405E 

06 

-97.215 

-26.98 

1992.32 

1 1*80000 

0.14082 

1.45 

2.971E 

09 

7 

2.337E 

06 

38. 162 

-8.54 

1994.46 

14.76667 

0.1430^6 

1.43 

3.259E 

09 

8 

2.805E 

06 

6.075 

-5.70 

1994.27 

1 6.89999 

0.12839 

1.41 

3.710E 

09 

9 

3.506E 

06 

-25.009 

-14.22 

1995.21 

18.96666 

C.l 1 709 

1.45 

4.749E 

09 

10 

3.504E 

06 

-56.094 

-22.73 

1996.85 

21.03331 

0.11561 

1.44 

4. 193E 

09 

11 

2.730E 

06 

-87.680 

-25.57 

1997.54 

23. 13332 

0.13194 

1.43 

2. 9666 

09 
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Output Table 2 - Electrons 
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